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ABSTRACT 


This paper documents the development of a non-para- 


metric statistics package for the TRS-&80O microcomputer. 
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emphasis on mon-parametric hypothesis testing. 


The pregramming lancuace is TRS-S0O Level II Disk 
Basic. The nackace is compatible with any disk based 
RS-80 Medel [/ ffi comnatible micrccompouter. With 


modification of the screen displav and the disk acces 


i 


commands, the package is tramsportable to micrcecomoputers 


produced by other manufacturers. 


fae ea eri Calbwamae!l vVSis capability imelemented on a 
|jFotatively inexpensive svstem orovides a useful! tool to the 
student er the trained analyst. 
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: INTRODUCTION 


Statistical Analysis of data is utilized as a decision 
aid in most every field. Fersonnel in technical fields 
have been exposed to and require the use of many 


statistical analysis prccedures during the course of their 


fu 


career. In a career in which statistical analysis Gf data 


is required , but not on a daily basis. the knowledae of 


é 


often retainsd where 


ps 
if 


the orcper statistical orocedure 
the mechanics of implementing the oroccedure reauires 


extensive review. A software packace which eliminates this 


J 


review and handles the mechanics of imolementing the 


precedure is most valuable. 


Tre cerscnmal computer has evolved to the degree that 
it 25 ncw well suited for the imolementatican of statistical 


analysis precedures fer small to mederately sized data 


p--d 


mM ow rr 
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ient package which models parametric 


distributions and implements parametric hynoethesis testin 


was developed by Mr. R. F. Isbell and was the subject of 


his thesis work at the Naval Postcraduate Schocl. The 


“ 





purpose of this thesis is to provide a software package 
which implements the more commonly used non-parametric 
procedures, accurately approximates their distributions, 
and orovides a convenient data base manacer useful for both 


this package and that developed by Mr. Isbell. 
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A. PACKAGE CONTENTS 


The package may be ltegically divided into three 
Nartitions; the data base manager, non-parametric 
distributions, and nen-parametric hypothesis tests. The 
data base management ocroegrams cermit data entry and storace 
to be in a columnar format cr ina tabular format. The 
orograms ocrovide a convenient file editing capability as 


well as a flexible file printing capability. 


The mere commenly used non-parametric distributions 
and their inverse distributions are aporoximated. The 


distributions are: 


hs Wilecoxson siqned rank distributicn 
a Mann-Whitney distribution 
es Smirnov distribution 
4. Kolmoqorov distribution 
iis filliefors modification to the KE-S test for 
Normal Distribution 
i1lilisefors modification to the K-S test for 
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“~Bponential Distrioution 


hi-Sauare distribution 
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The package contains orocedures for implementin 


FOllcwing hypothesis teste: 


Wilecoxsen signed rank test 
Mann-Whitney test 

Smirnov test 

. Test for Normal Distribution 

a. Mean amd variance unknown 

Ne Mean known, variance unknecwn 
cea Mean and variance known 

est for Exponential Distribution 
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ay Farameter unknown 
tH Parameter known 
oe Contincency Table Chi-Square test 


PE mle ROeOMPUTER SYSTEM 


The package was developed on an LNW Model I! 
microcomouter system with 46K bytes of memory and one sS-i/4 
inch floppy disk drive. The package 1s totally compatible 
with any TRS-8SO Medel I or Medel ITI disk based 
microcomputer system and is not dependent on any unique 
features of the advanced TRS-&SO Model [/IITTI operating 


systems. 


The oroqramming language 1s TRS-&O Disk Basic 
impiemented by Micrcsoft. Compatibility with all TRS-&0 
Medel i/III systems is maintained by using only TRSDCS 
compatible commands. The package 1s not directly 
transportable to non-TRS—-S8SO compatible comeuters; however. 
the Microsoft interoreter is cemmenly used by other 
mamufacturers. the disslay commands would have to be 
altered to be compatible with micrecemouters oroeduced by 


other manufacturers. 


With any procram develcoment on a microcomputer, the 
programmer must be cegnizant cf its speed. memory, and 


Seecision limitations. Te achieve maximum utilization of 


a, 


16 





the computer memory, the floppy disk 1s used as a virtual 
memory device pulling in routines on an as needed basis. 
This method permits a minimum of S090 data cbservations ta 


be resident in the computer memory. 


The speed of execution 1S a major concern when using 
an interpreted language. fe achieve an accentabies 
execution speed, asymptotic distributions were used in liesu 
Of exact distributions and faster but more memory intensive 


sort routines were routines were used. 


Precision is certainly a concern in this package; 
however, single precision variables proved to be adequate 
in most cases. Double precision variables were used only 
in the aoproximations to the Wilcoxson and the Mann-Whitney 


eretributions. 
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Ili. DISCUSSION Or FROGRAMS AND ALGORITHMS 


A. STRUCTURE 


The package consists of ten programs, 1.e. a control 
module, two data I/0 programs, a sSrogram calculating 
distributions and six programs performing hypothesis tests. 
The choice of the program to be executed is made througn 
the menu driven control module. The user 1s free to 
execute any orogram in the package through the control 
modul2® with all variables being initialized as a new 


srogram is executed. 


All data 1/0 1s conducted through the two data [/0 
procrams. The data 1s entered to the hypothesis test 
programs through a menu of established data files which 


were created by the data [/0 programs. 


B. PROGRAM 1, MENU/EAS 


MENU/BAS is the control medule from which all programs 
are executed. The program consists of the master menu and 
two sub-menus. The functional flow of the program is shown 


in figure i. 





BLOCK DIAGRAM OF MENU/BAS 


MASTER MENU 
Data Base Manager 


Probability Distributions 
Hypothesis Tests 

























DATA FORMAT MENU 


{. Coluaar Data 
Tabular Data 


Execute 
NPSTAT/BAS 


a a a a OD ee OD ee ee ee tee ee ee 


1. #ilcoxsen Sianed Ranks Test 
=. Mann-Shitney Test 
5. Sairnay Test 






4. Norgal Test Distribution 
s. Exponential Test Distribution 
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Execute Contincency Table fest 


DEMBR/SAS 

















Execute Execute 


SMIRNOV /345 


NUTLL/BAS 





Evecute 
CHIS8U/BAS 


ESURE  1:. 





mee FROoGRar 2. DEMGR/SAS 


DBMGR/SBAS orovides the data [/0 fer the hysothesis 
test programs. The program is menu driven and offers the 
following selections: 

1. Create New Data Files 

2. Edit Data Files 

3. Purge Data Files 
4. Frint Data Files 
a 


7 Return to Master Menu 


special purpose data file named EDITDATA/DIF 15 usec 


‘b> 


by this precram to contain the filenames of the data files 
erpeated. It is used as a data file directory and is called 
by the hyoothesis test ocrograms to display the tilenames cf 
all current data files and to permit a convenient menu 


selection of the data tiles. 


Meee §4« «XMEW FILE CREATION ROUTINE 


The new file creation routine prompts the user for the 
number of files to be created. The user is then sromoeted 
Poem the disk drive number, 1.e 0 - 3, to contain the files. 
The data files on the selected drive are dicepliayed to 
prevent any undesired duplicate file naming. The user is 
then promoted for the new filenames. if a previousiy used 
filename is chosen, the user is alerted that a duplicate 
filename has been used and the user is given the ertion of 


renaming the file or overwriting the existing file. 


bob 
3 





After the files are named, the input screen is 
isplayed. The screen will display four files (columns) 
and five cbservations (rows). The rows and columns may be 


scrolled to display other observations as needed. 


The input area is located in the top-left section of 
the screen. Derressing ENTER transfers the input area to 


the bracketed destination area. 


There are several special function keys which must be 


defined: 

ENTER = transfers contents of the input area 
to the destination area. If the ingsut 
area is Blank, the destination area is 
moved to the right. 

Sorrows _ moves the destination arsa 

Shift arrow - moves the destinaticn area to its 
most extreme positian 

D = deietes the current row 

I = inserts a blank row immediately before 
the current row 

CLEAR -- Feturns to menu 

=) — writes the files to disk and returns to 


the menu 


The data is written to disk in a sequential file with 


the first field containing the number cf observations and 


ui 


eacn field thereafter containing the observations. Th 


oa 
0 
< 
(D 
of 
(0 
“| 


routine was desianed for maximum user convenience; 


p> 
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these convenience features will permit undesired "zero 
valued observations if the user is mot aware of the design 
of the program. Some precautions and guidelines to follow 


ar@s« 


1. ALI disolayed files will be assigned a commen 
value for the number of ocbservaticns. Therefore. only 
files which will have the same number of observations 
should be entered simutaneously. 


e . It is sometimes tempting to "play” with the arrow kevs 
to move the destination area around the screen. The 
number of observations is determined by the maximum row 
reached by the destination area; therefore. "playing" 
with the arrow keys may result in an unwanted value tor 
the number of observations. 


See |Digg EXISTING FILES ROUTINE 


+ 
ai 
(D 
“h 
pa 
wd 
(U 


editing routine operates similiarly te the 
new fils creation routine with only minor variation. The 
user is prompted for the disk drive containing the files 
and then prompted for the selection cf the files to be 
edited. The files are displayed on the screen and editing 


1S performed exactiy as described in the new file creati 


é 


Q 
a 


Atl 
(D 
a 
p- 
Ul 


memgeree. Prieocr to writing the files to disk, the us 


Given the oppertunity to rename the file. 


Sees FURGE FILES 


The data files are dicnplayed one by one and the user 


15 asked if 1t iss desired to erase the file from the disk. 





All erasing of data files should be performed from the 
purge routine and not from the operating system routines as 
the purge routine also removes the filename from the data 


file directory, EDITDATA/DIR. 


C.4. PRINT DATA FILES 


This routine permits fermatting the data files for 
orinting. The user is promoted for the files to he 


oOrinted, the printer line width, column headers. amd civen 


‘rt 


he sotion of numbering the obtervatisns. The maximum 
number ef spaces between columns is calculated and 
displayed. The user 1s prompted for the desired number of 
snoaces between columns. Finally the user is asked for the 


report title and the 1198S are printed. 


D. FPRCGRAM 3, CTABLE/BAS 


CTABLE/BAS is used te create tables used in the 
contingency table test orccram. The oproaram operates 
Similiarly to the DEMGR/BAS orcoram with the differences 
listed below: 


a. The data file names are stored in CTABLE/DIR versus 
EVITDATA/DIR. 


ee In the new file creation routine, the user is 
Promoted for the number cof rows and columne in the 


table. 


s- The "DBD" and "I" keys have no function. 


f= 
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E. FROGRAM 4, WILCOX/BAS 


WILCOX/BAS performs the Wilcoxson Signed Rank test 
using data stored in data files created by DEMGR/BAS. The 
user 15 prompted fer the two data files tor the test. The 
test statistic 1s computed and the user is orompted for the 
hypothesis testing problem, 1.2. cne sided or two sided. 
jhe test statistic, sample size and p-value are displayed 
and the user 1s prompted for the significance level of the 


test 


— 
é 
‘ 


. he significance level is compared with the o-value 


and acceptance cr rejection of the null Rynpothesis 15 


determined. 
1. Wit COXSON SIGNED RANK TEST 


The data consists of two paired samples of n 
ebservations, ( X4s¥q)q (Mas Yo) sees f 
The difference d:= x:-y; 1s computed fer each cf 
the n pairs. Pairs with a difference of zero are omitted. 
Bee absoittte differences, !d,;!, of the remaining sairs 


are ranked with rank 1 assigned to the pair (x;.,.y;) 


with the smallest absolute difference. If one or more 
Pairs have the same absolute difference. the average rank 
1S assigned to the pairs. The following test statistics 
are computed: 


+ ‘ : 
T° = Sum of ranks with sositive Gs 


3 
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VY 


in 


T = Sum of ranks with negative oc 


pond 
CO 





Te eranich 7) 


The assumptions for the test are: 
ieee tine GiStributicon of each d; is symmetric 


Senie pairs (4 ,5¥,7 cCenstitute a bivariate random 


sample. 


The hysothesis may be stated in ome cof three ways: 


sets eee ee Ges 
Hy, 2 EtX) x ECY) 
Bo Ho = EtX) 2 EtY) 
Bien ss Et Y) 
C. Hy, : EfX) 2 ECY) 
re ee ee Et yD 


The decision rule for the three hypotheses is to accent 
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See Fitest statistic « 73 * Significance Levei/ss 
Bee (feet Statistic =< 1) Significance Level 
Pei teet Statistic =< fF) SeeSionificance Level 


Pees JDISTRIBUTICN THEORY AND ALGURITH! 


The distribution of the Wilcoxson statistic. T. may be 
found by randcemizing the signed ranks. Since there are two 
pessible signs, + and -, andn ranks, the number cf 


combinations in the randomization is 2. Thee ae eC 


SeGemnOoinaticons and counting those with the sum of ranks 
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less than or equal to the test statistic and dividing By 
2?) Such a method is cbhviously siow for all but very 


smali sample sizes. 


A number of aporoximaticns to the exact distribution 
have been proposed. The most accurate of the proposed 
aporoximations is the Edgeworth arporoximation which is 
documented in References 1. 2, and 3. This approximation 
is based upon the fact that the distribution of Tis 
asymototicaily nermal with mean and variance: 


Ur = mint+i)/4 


qn* = nintt) (Pn4+1i)/24 
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normalized Wilcexson Statistic, z, is defined: 
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ne Edgeworth exnansion to terms cf size i/vn~ is 
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wnere Fiz) is the cumulative standard nermal distribution 
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F.2. DISTRIBUTION THECRY AND ALGORITHM 


The distribution of the Mann-Whitney statistic may be 
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mee) 6 FPROGRAM 7. NMLILL/BAS 


MLILL/BAS performs a test for normal distributicn ona 
data set stored in a data file created by DEMNGR/BAS. The 


test may be performed for the three cases 


ae mean and variance unknown 


ae mean known and variance unknown 


Be mean and variance known 
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The user is proemoted fer the data tile to be tests 
and prompted fer the specific case to be tested. In the 
cases where a parameter is known, the user is ocromoted far 
the known parameter. The test statistic. samrele size and 
D-value are displayed. The user i¢c prompted for the 
Significance level of the test. The significance level is 
compared with the o-value and acceptance or reijectisn of 


the null hynothesis is determined. 
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The assumption fer the test is that the sample is a 
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For Case 3, where the mean and variance are completely 
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sample size and p-value are displayed. fhe user is 
orompted for the significance level of the test. the 
significance level is compared with the p-value and 
acceptance or rejection of the null hypothesis is 
determined. 
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the distribution was approximated using Monte Carle 
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considered to be satisfactory if no Ei is less than i 
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The algorithm used in this package for the Chi-Square 
distribution was extracted from the thesis of Mr. R.F. 
Isbeli. The details of his approximation to the 
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Pefee., a Correction factor sronosed by F. Yates in Reference 
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MPSTAT/BAS contains seven non-parametric probability 
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cay OO 


y 


pei actos, 

pe oi oO 
aS - 5454 
22044 2SOSo 
eo? 27 oS 
- 2580 c2oe4 
ee 7. mea ay 
2474 me 
ee Sy 


@ ae So” 


1 
> 
4 
Ss 
oO 
7 
8 
9 
Q 


“am / aA 
Se se 2ae 


CODON & OB OH hoe 


foo59) 
so7S 
eee a 
ao A 
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TABLE 4 


KOLMOGOROV DISTRIBUTION ACCURACY COMPARISON 
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TABLE 5S 


Seer Srenomieinmel ity (251 DISTRIBUTION 
ACCURACY COMPARISON 


= Eywact 








TABLE 65 


Pee Pe oee Ae shemMTieiity FES! DISTRIBUTION 
ACCURACY COMPGRISUN 


Sampie 
Size 





Ci 
( 
8 
o 
( 


Re) al ee io so ite mee te? 
oi Ne ee eo BSS a ioe 
as ~ i447 FLOSS eee Gee oe 


pot 
i 
tod 
pense 


pub 
Cr- 
an) 
bob 
p= 
Ci 

ch] 
id 


>—L 
ofa 
fa 
| 
po 
(4 
(0 
(0 
p> 
Ci 
(3 
Q 


f+ 
i 
“4 
i 
>p~> 
Ci 
C3 
Q 
-~ 
cn 
ioe 
cy 


aS foot eee 


a 
p~> 
( 
U 
p= 
a 
pb 
tJ 
a c is ® 
Cs 
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bd 
nd 
a 


Si) oc os 
et 0914 SOLS ele Se BS ae es 
150 0862 eS ye oS Se 








PROGRARP 2 = MENU, RAS 
10 CLEAR 350 
20 DEFINT I-K : DEFSTR Z 


30 LBS=CHRS (123) 
BL$= 
40 CLS:PRINT TAB(10) 


RBS=CHRS (125) MBS=CHRS (176) 
Cain? (252) 


"NON-PARAMETRIC STATISTICAL ANALYSIS 


SOFTWARE" 


30 
&0 
79 
of 
=t@) 
9O 
100 
110 
120 
130 
140 
150 
160 
170 
180 
199 
=O00 
Z19 
220 
20 
240 
eet) 
3O0G 
319 
520 
53G 
340 
NOT 
DEP, 
a0) 
19000 
1010 
1020 
1930 
1040 
1050 
1040 
1970 
1G80 
1100 
11190 
Bee 
1150 
1140 
1150 


o 
£ 
oS 


PRINT 
PRINT 


PRINT 


TAB(S1) “by" 

TAB(21) “Robert Lee Zangmeister" 

PRINT TAB(C1S) “Submitted in partial fulfillment 
PRINT TAB ‘17) "requirements for the decree of" 
PRINT TAB(19) "MASTER OF APPLIED SCIENCE” 

PRINT TAB(29) “from the" 

PRINT PRINT TAB‘19) "NAVAL POSTGRADUATE SCHCOL" 
PRINT TAB(C25) “December 1982" 

PRINT:PRINT TAB(C17) "Press SPACE BAR to continue" 
ZI=INKEYS feet THE 150 

ELS PRINT TAB(2Z0);"PROGRAM SELECTION MENU" 

PRINT 23128," 1) Data Base Manager" 

Prim? "(s) Probability Distributions" 

PRINT “(3) Hypothesis Tests" 

KP=996 mS=6 GOSUB S00 

ON I0 GOTO 200,250, 1000 

GESSeRIMNT “DATA BASE OPTION’ :PRINT 

PRINT "¢1) Column Format Data Entry" 

minim 6 C2) Table Format Data Entry” 
RP=896:KS=2:GOSUB S00 

IF i090 1 THEN RUN "DBEMGR/BAS" 
RUN “NFPSTAT/BAS" 

REM * & * KEYBOARD ENTRY ROUTINE * 

aa wer tS, OPTIGNe DESIRED “:LBSsibSs RES: 

ZI=INKEYS IF Zi="" THEN S20ELSE TO=VAL (ZTE) 

PRINT oKP+21, Zis: FOR I=1 TO 194 NEXT 

Peeve. OR LSekS THEN PRINT OEP, BLSs: PRINT oP, "Se 


RINT 
he" 
RINT 


Bee RUN)  CaeaSley BAS* 


A VALID OFTICQN keees: FOR f=1 7O 1900 : NEXT : PRINT 
mete : GUIG Silo 

RETURN 

Seem tebe eO) s YP OTHESIS TEST MENU": PRINT 

PRINT 3128,"41) Wilcexson Siqned Ranks Test" 

PRINT "(2) Mann-Whitney Test" 

PRINT "{3) Chi-Square Contingency Table Test" 

PRINT "(4) Smirnov Test" 

PRINT "(35) Test for Normality” 

PRINT "{6) Test for Exponentiality” 


KP=8946:KS=6:GOSUB SG0 

eee O GETS 11Oov1iio, 12120,1150,1140,1150 
RUN "WILCOX/BAS" 

RUN “MANN/ BAS" 

RUN “CHISQU/S3AS" 

RUN “SMIRNGV/BAS" 

RUN “NLILL/BAS” 

RUN "“ELIL!E/BAS" 
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rF-RoGRArM ~~ = DEMGR 4 RAS 


10 CLEAR 1000 

20 ON ERROR GOTO 19000 

50 DEFSTR Z 

me DIM FPSd{9S),FiS (il) ,MSc1l> 

50 LES=CHRS (1235) : RBS=CHRS (125) : MBS=CHRS (176) : BLS=CHRS (252) 
60 CLS 

70 PRINT TAB(20)3;"DATA FILE MANAGEMENT PROGRAM" 

BO PRINT: PRINT: PRINT "MENU: " 

90 PRINT "(1) Create New Data Files" 

100 PRINT "(3) Edit Existing Files" 

110 PRINT "(S3) Purge Existing Files" 

120 PRINT "¢4) Print Data Files" 

130 PRINT "{3S) Return to Master Menu" 

140 KP=896:KS=5:GOSUB 180 

150 Pee Seotetneeol., 51, S50, 540,170 

1460 GOTQ 40 

170 RUN "“MENU/BAS" 

180 FEM * x* KEYBOARD INPUT SUBROUTINE 

mo PRUNT aEP+S, "OPTION DESIRED ":;LBS;MBS;:RES; 

feo ZI=INIEYS: IF Zi="" THEN 200 ELSE IG=VAL (ZI) 

oem ONT agKP+t21,Z1s:FOR f=1 TO 100: NEXT 

220 IF IO<1 OR [Os2KS THEN PRINT OKP,BLS3:FRINT OEP," *% 
NOT A VALID OPTION x*x"s: FOR T=1 TQ 1000: NEXT: PRINT 
fee, Siw sGOTO 1979 

250 RETURN 

240 DNS=INKEYS$: IF DNS=""" THEN 240 ELSE I=VAL (DNS) 

230 IF I20 OR Is3 THEN 240 ELSE RETURN 

=50 *SUB TO PERFORM NEW DATA FILE CREATION RCUTINES 
=7Q09 CLS : INPUT "How Many Variables are to be Used 
Seti) "se: CT=NE: RT=1 

280 IF NF<i OR NFSitl THEN 270 

290 GOSUB 14000 : GOQSUB 12000 : N=1 : CR=1 : GOSUB 6020 
S00 CLS:GOSUB 15050: RUN 

510 GOCSUB 350 : CR=1 : GOSUB 6020 

529 GOSUB 12000 : RUN 

eee COSUB SOGO : RUN 

340 GOSUB 3bO : GOSUB 4000 : RUN 

300 REM x * * READ MULTIPLE ARRAYS x *X 

360 K=0 : DLs="/" 

of/9 CLS:PRINT "Disk Drive Containing Files (0-3)7" 

580 GOSUB 2490 

S979 SL=1715:GOSUB 1200: IF SL=0O THEN GOTO 370 

4090 GOSUB 1400 : FRINT "" 

410 PRINT:IF I0 = 2 PRINT "Enter the Fiies To Be Edited 
fe.g. 1/3/4)"3 ELSE PRINT “Enter the Files To Be Printed"; 
429 INPUT BS :Pi=0 

450 P2=INSTR(1,58%, DLS) 
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440 
450 
440 


IF P2=0 THEN NV=VAL (BS) : 
NV=VAL (LEFTS (8S, P2-1)) 


IF NV>FO THEN PRINT: PRINT" 


R=) GOTO 460 
BS=RIGHTS (8S,LEN(BS)—-P2) 
* * XK INVALID FILE # xk 


*":GOTO 410 


470 
420 
490 
300 
ANOot 
310 
Zor = (s 
320 
Soo) 
S340 
3UO0 
360 
370 
380 
39 
6ao 
610 
620 
630 
540 
5350 
660 
1900 
1019 
Lo2D 
1030 
1040 
1050 
1060 
10709 
1190 
11190 
1120 
1130 
1200 
1210 
ee) 
1230 
1240 
1250 
1260 
1270 
1280 
1290 
1300 
1310 
1400 


K=K+1 IF K=12 THEN NF=1ii 

Fid(K) =F S$ (NV) 

IF Pi=0 THEN 430 
BS$=""SPRINT:PRINT "Do You Wish To Read Files From 
her Disk (Y/N)? 
ZI=INKEYS: IF ZI="" 
OTO 370 
NF=KsCTH=K:2:PRINT 
MM=-—1 
FOR I=1 
MisceNe st  ,1,Fis (I) s 
IF M > MM THEN MM=M 
NEXT I 
RT=MM 
GOSUB 12000 
FOR J=i TO NF 

PRINT "READING IN FILE 

OFEMe “we”. 1. Fisid) 

FOR I=1 TON 

CLGSE 1! 

NEXT J 
N=MM =: RETURN 
REM DIRECTORY CHECK 
IF FO=0 THEN RETURN 
P2=0 
FOR Pi=1 TO FO 
Sens toe i,FSetPi-P2) 
Geese 
NEXT Pi 
IF P2=0 THEN RETURN ELSE GOSUB 
REM * * * CORRECT DIRECTGRY 
IF Pi-P2=FO THEN 1130 
FOR I=Pi-P2 TO FoO-i FPS(IT)=F$(1+1) 
FO=FO—1 PD=P2+1 RETURN 
REM * * * READ DIRECTORY 
CLS 
FS="EDITDATA/DIR: "+DNS 
SPEN “H",1,F% 
INPUT #1,F0 
IF FO=0 THEN 1299 
FOR I=1 TO FO 
INPUT #1,F$¢(1) 
NEXT I 
Beg Se 
SeGSe 1 
RETURN 
REM «x x & 


Cane las 0 


as 
Sereno lo else slr Zr y¥" PRINT 
nN 


10 -NE 


INPUT #1, GEeSE 1 


PRINT 

—-> "sF1$(J) 
INPUT #1,N 
INPUT #1,X(I.J) 


NEXT I 


x xX 


2400 
x 


RETURN 
x 


NEXT I 


x x 


60SUB 19000 


DISPLAY DIRECTORY x £ 
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1419 
1420 
1430 
1440 
1450 
1460 


PRINT 


1470 


GES 

IF FOQ=0 

FOR I=1 
FOR 


THEN 1510 

TOSEO STERA4 

J=1 TO 4 

Teei+td-1i7PO THEN 1510 

IF J=1 AND LEN(FS(I+3-1))<¢12 AND I>1 THEN 


PRINT TABCITGSRI-1))35 °<("SSIRS¢(I+3-1)3") 


ere (I+J—1); 


1480 
1490 
i500 
1510 
2500 
2505 
aoe 
240 
23545 
2360 
2370 
S400 
2410 
3420 
3439 
3446 
S450 
35460 
3470 
2480 
SOOO 
So19 
SILO 
Sa 

3Su4O 
SOO 
35600 
5610 
S620 


NEXT J 
NEXT I 
PRINT 
RETURN 
REM x * & 
CLOSE 1 
F$="EDITDATA/DIR: "+DNS 
OPEN "0",1,F% 

PRINT #1,0 

CLOSE 1 

RETURN 

REM «* * * ENTER NEW FILES «* 
F$="EDITDATA/DIR: "+DNS 

OPEN "0",1,F% 

PRINT #1,F0 

FOR I=1 TO FO 

PRINT #1,F8(I)3","3 
NEXT I 
CLOSE 1 
RETURN 
REM «x * X 
EX=0 
FOR I=1 TO FO 

IF FS=F$ (I) 
NEXT I 
RETURN 
REM «x * * 
CLS 
INPUT 


CREATE DIRECTORY x & 


FILE NAME CHECK x x 


THEN EX=1 : RETURN 


DIRECTORY DISPLAY * X 


"DISK DRIVE NUMBER (0,1,2,2) 


We oP UND DIRECTORY 


OF"; DNS 


3630 


3640 


IF VAL (CDNS)<0 OR 
SUSUe L200 5 


VAL (DNS) +3 THEN 
GOSUB 1400 : 


mo.) 


PRINT "" : INPUT "PRESS 


ENTER TO CONTINUE"; Ji 


2650 
4000 
4005 
4010 
4920 
4030 
4040 


RETURN 
REM * * * PRINT DATA FILES * ¥ 
CLS : INPUT "Input Printer Width";FuW 
PF=0 : ML=-1 : DL$="'!" : TL=-1 
FOR I=1 TO N 
FOR J=1 TO NF 


Pao itesine «A Ces J) )) 2 ME THEN 


ML=LEN(STRS(X¢(I,9))) 


4950 


NEXT J 
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4060 
4061 
4070 


NEXT If 
ML=ML+1 
GES: PRINTSPRINT 


"Enter Column Headers. Separating 


lines with" 


#080 PRINT"the delimiter *!'* to produce a two column 
header" 

4100 FOR {I=1 TO NF 

4110 PRINT:PRINT "Enter the Column Header for Column"; ! 
4120 INPUT MS$(T) 

4130 NEXT I 

4140 MS=-1 

4150 FOR [I=1 TO NF 

4160 P1=0 

4170 P2=INSTR(P1+1,M$(1I), DLS) 

4190 =P 2-e1— 1 

4190 IF P2=0 THEN L=LEN(MS(1T))—-F i 
4200 ITF L>MS THEN MS=L 

4210 IF P2=0 THEN 424090 

4220 Pi=P2 

4230 GOTO 4170 

4240 NEXT I 

4250 IF ML>MS THEN MX=ML ELSE MX=MS 
4260 FOR T=1 TO NF 

4276 P1=0 =: NL=0 : Ts="" 

4299 P2=INSTR(P1+1,M$(1I).,DLS) 

4300 L=P2-P1i-1 

43190 IF P2=0 THEN L=LEN(M$(1T))—-P 1 
4320 SF=INT(CCMX-L)/2) =: SB=MX—-L-SF 
4530 TS=TS$+STRINGS (SF, " 


"Y+MIDS(MS(I) ,Pit+i1,L)+STRINGS (SB, " 


hee 


#555 NL=NL+1 = IF NLSTL THEN TL=NL 

4346 IF P2=0 THEN M$¢(T)=TS =: GOTO 43570 

4350 Pi=P2 

4360 GOTO 4290 

4370 NEXT I 

4380 FRINT:PRINT "Would you like to number the observations 
masN)? "s 


43585 AS=INKEYS:IF Ags="" 
4390 IF Aag="y" 


THEN 4582 
THES? Pi=S ELSE PI=0 


ELSE PRINT AS 


4400 MS=INT ( (PW-MXS¥NF-PI) /NF) 


4410 PRINT:PRINT"The maximum # of 


spaces that can be used 


between columns is";MS;"." 
4420 PRINT: INPUT"Enter the number cf spaces between 
columns"; SP 


4445 IF SP<O THEN INFUT 
N"sTS 
4450 
4455 
4472 
4474 


"TO RENAME COLUMNS TYFE Y ELSE 
>: IF F#="Y¥" THEN 4010 ELSE RETURN 

IF SP<0 OR SPOMS THEN RETURN 

PRINT:INPUT “Enter Report Title";:TSs 

Pi=(NF-1) * (MX+SF)+PI+MX—-1 


meine (RP l/ 2) 5 PL=PW-P1 =: IF LEN(TS) >PW THEN 


TS=LEFTS (TS, FW) 
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4476 IF LEN(TS$) -P2 


TREN Pi=INT ((BW-LEN(TS))/2) 


4480 LPRINT TAB(P1);TS 


4490 LPRINT TABC(P1i); : FOR Y=1 TO LEN(TS) =: LFPRINT 
Pee sexy f =: LPRINT 

4500 LPRINT "" : LPRINT "" 

4505 FOR J=i TO TL 

4310 FOR T=i TO NF 

4520 LPRINT 

pee (I-17) 2 (MX+SP) +P1I) sMIDS CMS C(I), CJ—-1)X+1,MX); 
4530 NEXT I 

4340 LPRINT "" 

435350 NEXT J 

4360 AS=""": FOR T=1 TO MX =: AS=AS+"—" = NEXT I 
4370 FOR I=1 TO NF 

4580 LPRINT TABC(CI-1) (MX+SP)+PI) 3s AS; 

43590 NEXT I 

4600 LPRINT "" : K=PI+MX—-1-INT ((MX—-ML) /2) 

4610 FOR I=i TON 

4620 foe —o SEEN DEPRINT STRS¢1)-": "*s 

44550 FOR J=1i TO NF 

4540 LPRINT 

em. (J —1) SXMATSP) —LEN(STRS(X(CI J) DD +K) SS XCIST) S 
4620 NEXT J 

4560 LPRINT "" 

4570 NEXT I 

4580 RETURN 

2000 REM * * * PURGE DATA FILES x & 

go10 CLS 

39020 PRINT “Disk Drive Containing Files to Be Purged 
me—35)?" 

3050 GOSUB 240 

9040 SL=S5010 : GOSUB 1200 : IF SL=0 THEN GOTO S010 
2045 K=0 

3046 PRINT TAB(20)3" kXe*4K PURGE DATA FILES **xx" 


2030 
3060 
JO61 
sO62 
3064 
sO90 
Se75D 
So79 
OO 
=1ii0 
= 20 
e220 
3140 
5000 
6003 
6010 
601i 


FOR I=1 TO FO 


AS=H=""SPRINT FSCI-K) 3" Poetieeiae CY /MN7E) "> 
ZI=INKEYS: IF ZI="" THEN S061 

ieee — VO” THEN PRINT ZisGOTO 3095 

ieee —  Y PREM PRINT Z2ZisGOSUB 3Sigo0 ELSE PRINT “N" 
NEXT I 

GOSUB 3400 

RETURN 

KILL FSCI-K) 


Beeet—FO THEN 3130 
mem J—i-k TO FO-i = FS(J)=FS(3J+1) = NEXT Jd: 
FO=FO-1 3: K=K+1 
RETURN 
REM **%x EXAMINE DATA x 
INPUT CR 
IF CR < 1 THEN CR=1i 
IF CR > N THEN CR=N 


a9 





6020 ES=CHRS(31) =: NS="" = P1=0: CLS 

6025 Ci=i:Ri=i:Bi=1 

6020 CR=R1i:CC=Ci 

6035 GCSUB 9500 

6040 GOSUB 6300 : GOSUB 6350 

6045, P1=0 

6050 ASB=INEKEYS : IF ASG=""" THEN 6050 

60SS5 IF AG=CHRS(91) AND Pi=0 THEN NS="": GOSUB 6500 : GOTO 
6050 

6060 IF ASG=CHRS(27) AND P1i=0 THEN Ns$="": CR=1:GOSUB 
46520:GOTO 6050 

6065 IF AS=CHRS(10) AND P1i=0O THEN N#@=""; GOSUB 7000 : GOTO 
6050 

6670 IF AS=CHRE(25) AND P1=0 THEN N$="": CR=RT: GOSUB 

7020: GOTO 60350 

6075 IF AS=CHRE(3B) AND Pi=0O THEN N$="": GOSUB 7500 : GOTO 
60590 

6080 IF AS=CHRE(9) AND P1i=0Q THEN NG=""; GOSUB 8000 : GOTO 
60350 

6085 IF AS=CHRS(24) AND Pi=0 THEN N$="":CC=1:GOSUB 

7350: GOTO 6a50 

5090 IF AS=CHRS(25) AND FP1=0 THEN NS="": CC=CT: GOSUB 
8030:GOTO 6050 

6095 IF (AS>"/" AND AS<":") CR At="_." OR ASH="—-" OR AG="+" 
OR AS="E" THEN Pi=1i : GOSUB 9000 : GOTO 6050 

6100 IF AsS=CHR#(13S) THEN Pi=0: GOSUB 8500 : GOTO 60350 
6105 IF At=CHRE(S8) AND Pi=1i THEN GOSUB 90290 : GOTO 6050 
6110 IF A¢="D" AND Pi=0 THEN NS$="":;GOSUB 6320: GOSUB 
19600: GOTO 6950 

6115 IF 4$="I" AND Pi=0 THEN NS="":=GOSUR 6320: GOSUB 
190500:GOTO 6050 

6120 [F AS=CHRS(S1) AND P1=0 THEN RUN 

6125 IF AS="9" AND P1=0 THEN N=RT : RETURN 

61350 GOTO 6050 

6300 REM * * * CLEAR ENTRY AREA * X 

6310 CC=1 

6320 PRINT 22," ie 

65350 PRINT » 2, CHRS (9S) 

6540 RETURN 

6350 REM * * K BRACKET ENTRY & 3 

6360 K=324+128% (CR-R1)+15% (CC-Cl1) 

E370 TY=15260tK 

Seo FORE TY,185 : POKE TY+t14,187 

6390 RETURN 

6400 PORE TY,128 : POKE TY+14,128 

6410 RETURN 

6300 REM * * * PROCESS UP ARROW x x 

6310 CR=CR-1:I1F CR<1 THEN CR=1:GOTO 6540 

6329 IF CR<RI THEN Ri=CR: RL=R1+4:6GO0SUB 2500 

6550 GOSUB 6400:GOSUB 6350 


6540 


RETURN 





7900 
7910 
7020 
7030 
7040 
7030 
73500 
7210 
Fuso 
7350 
7I40 
7330 
Boon 
8010 
8620 
8650 
3040 
8050 
8500 
8510 
8520 
a530 
8540 
8550 
8560 
8570 
8280 
8590 
84600 
S610 
FQa00 
$010 
CHRS 
P0209 


9OSO 
S220 
9210 
9220 
$230 
9240 
es () 
92460 
9270 
9500 
oo 1a 


eo) : 


REM * * * PROCESS DOWN ARROW «* x 
CR=CR+1 

IF CR>RL THEN RL=CR:R1=RL—-4:GOSUB 6400:GOSUB 
IF CR > RT THEN RT=RT+1 

GOSUB 6400:GOSUB 6350 

RETURN 

REM * * & 
CC=CC-1 

IF CC<1 THEN CC=1:GOTO 7550 

IF CC<C1 THEN C1=CC:CL=C1+3:GOSUB 9500 
GOSUB 6400:GOSUB 6350 

RETURN 

REM x & * 
CC=CC+1 

IF CC > CT THEN GOSUB 9200 

IF CC>CL THEN CL=CC:Ci=CL-3:S0SUB 9500 
GOSUB 6400:GOSUB 6350 

RETURN 

REM «* & x 
GOSUB 4320 
IF N$="" THEN GOSUB 8000 : 
IF CR>RT THEN RT=CR 

¥ (CR, CC) =VAL (NS) 

N= us 94 

=325+128% (CR-R1)+15* (CC-C1) 
PRINT OK+1,X(CR,CC); 

RETURN 

K=K-2:PRINT 3K,CHR$ (128); 
PRINT USING 2Z2;X(CR,CC); 
RETURN 
REM * * & 
NS=NS+AS : 


93500 


PROCESS LEFT ARROW xX x 


PROCESS RIGHT ARROW x & 


CHANGE VALUE IN FIELD * & 


RETURN 


FROCESS NUMBER ENTRY 
PRINT a2, “ ee 
RETURN 
IF LEN(NS) 20 THEN NS=LEFTS(NS.LEN(NS)—-1) : 

"s PRINT D9 2,NS 
RETURN 
REM «*« * xX 
Ce=1 : Ci = 
CR=CR+1 
IF CR>RE THEN RL=CR : 
GOSUB 6400 

IF CY > 4 QR CR = RL 
IF CRERT THEN RT=RT+1 
RETURN 
REM * * * SCREEN DISPLAY * xX 
Pree gies. CHRSt las) = IF CT<Cit+t4 THEN CL=CT ELSE 


x 


PRINT 2 2,N$::FRINT 


PRINT 3 2," 


LAST COLUMN PROCESSING « 3 


Cee eS 
R1i=RL-4 


feet "6GsuE FSeo 


ee—C1+S 
9320 RL=R1+4 


ease FOR I=Ci-1i 


9540 


Teoer— 1 
aii 197+ i Geet l+i-Ci), STRINGS (6, CHRS (143) )-:2PRINT 


fe) 20) SPRINT CHRS(65+1)s:PRINT CHRS (128) 5:PRINT 


a 





STRINGS (6, CHR$(143)) 

9550 NEXT I 

9560 PRINT 3195,CHRS(170):PRINT 2259, CHRS(170) 

9570 FOR I=Ri-1 TO RL-1 

9580 IF I > 8 THEN K1i=0 ELSE Ki=1 

9590 PRINT 93270+128% (1+1-R1)+K1,STR$(I+1)3:PRINT CHRS$(179) 
9400 FOR J=Ci-1 TO CL-1 

9610 PRINT 2526+15*(J+1-C1) +1294 (I+i1-Ri),X(I+1,J+1) 

9420 NEXT J 

94350 PRINT 2387+128*(I+1-R1) ,CHR$(170) 

9640 NEXT I 

9650 RETURN 

19000 REM * x * DELETE ALINE * x 

10010 FOR I=CR TO RT-1:FOR J=1 TO CT:X(I.J)=X(I+i,J) =: NEXT 
J : NEXT I 

10020 FOR J=1 TO CT: X(RT,J)=0 : NEXT J 

10030 RT=RT-1 : IF Ri > 1 AND RL > RT THEN R1i=Ri-1 : 
CR=CR-1 

190040 GOSUB 9500 : GOSUB 46250 

19050 RETURN 

10500 REM * * * INSERT A LINE «x x 

16510 FOR I=RT TO CT STEP -1 : FOR J=1 TO CT : 
X(I+1,J)=X(I,J) +: NEXT J : NEXT I 

10520 FOR J=1 TO CT : X(CR,J)=0 : NEXT J 

10530 RT=RT+1 : GOSUB 9500 : GOSUB 4256 

10540 RETURN 

11000 REM «x x * VERIFY FILENAME «x x 

11010 FW=0 

11020 P1=INSTR(1,S$,"/") 

11050 IF Pi <> O THEN IF LEN(S%$)-Pi>3 THEN FW=1 : RETURN 
11040 IF Pi >8 THEN FW=1 : RETURN 

11050 IF Pi=0 THEN IF LEN(S$)?@ THEN FW=1 : RETURN 

11060 IF Pl <> O THEN IF ASC (MIDS (S$, F1+1,1))<45 OR 

ASC (MID$(S$,P1+1,1))2>90 THEN FW=1 : RETURN 

11070 IF ASC(LEFTS(S$.1))<465 OR ASC(LEFTS(S$,1))>90 THEN 
FW=1 : RETURN 

11080 IF Fi=0 THEN P2=LEN(S$) ELSE P2=P1-1 


11090 P==1 

11100 FOR I=PS TO P2 

11110 AS=MID$(SS,I,1) 

120 IF (ASC(AS)>47 AND ASC (AS) <58) OR (ASCLAS) 264 
AND ASC (AS)<91 ) THEN 11140 

SO FW=1 2 RETURN 


11140 NEXT I 

feo 1- PS=i AND Pil <> 0 THEN PS=FPi+1 s P2=LEN(53) =: GOTO 
PerOd 

11160 RETURN 

12000 REM * * * DIMENSION ARRAYS X* % 

120160 NV=MEM-S00O : NV=INTONV/ (4% (NF+1)))-1 

feos? IF MM > NV THEN PRINT “INSUFFICIENT SPACE TO READ IN. 
"sINPUT "PRESS ENTER TO CONTINUE"; AS 


Ch 
tJ 





12030 DIM X(NV,NF) 

12040 MN=NV 

120350 RETURN 

135000 REM x * X RENAME FILES ROUTINE «x x 

mooi CLS 

13020 PRINT "Disk Drive Te Contain Files (0-3)°?” 

13030 GOSUB 240 

135040 SL=8240:GOSUB 1200 : GOSUB 1400 : FRINT "" 

13050 FOR J=1 TO NF 

13060 PRINT "Variable #"3J; "Currently Uses Filename 
fae a 1S) 

13070 INPUT “New Name";S%3 : GOSUB 11000 

12080 IF FW=1 THEN PRINT "INVALID FILENAME" : GOTO 
13060 

13090 FS=S$ =: GOCSUB S500 

13190 IF EX=1 THEN PRINT:FRINT "A File Currently Uses 
the New Filename.":PRINT "Do You Want to Use it Anyway 

er Noe" ELSE GOTO 13120 

13110 AS=INKEYS: IF AS=""" THEN 13110 ELSE PRINT AS: IF ASKS 
oy )| 6THENM: 13060 

foi 20 Fisi(J)=FS 

13130 NEXT J 

12140 RETURN 

14000 REM x * &* NAME FILES ROUTINE x X 

14019 CLS 

14020 PRINT "Disk Drive To Contain Files (0-3)?":60SUB 240 
14030 SL=8330 : GOSUB 1200 : GOSUB 1400 : PRINT "" 

14040 PRINT 

14050 FOR J=1 TQ NF 

14060 PRINT "FILENAME FOR VARIABLE #"3;J39 : INPUT SS : GOSUEB 
11000 

14070 IF FW=1 THEN PRINT "INVALID FILENAME" =: GOTO 14060 
14080 F$=S$ : GOSUB 3500 

14090 IF EX=1 THEN PRINT “FILENAME IS ALREADY USED. TYPE 
Y TO USE NAME ANYWAY ELSE N"s: INPUT AS : IF AS <> "YY" THEN 
14060 

14160 Fiti(J)=FSs 

14110 NEXT J 

14120 RETURN 

13000 REM x kK k* EXITING ROUTINE FOR OPTIONS 1 & 32 k x 
135010 CLS 

15020 PRINT "Do You Wish To Rename Files (Y/N)7?" 

15020 ZI=INKEYS: IF Zit=""" GOTO 15030 

15040 IF ZI="Y" THEN GOSUB 13000 

Bees SL=eso0 - GOSUB 1200 : GOSUE 1406 : PRINT "" 

15060 PRINT: PRINT 

13070 FOR kK=1 TO NF 

15080 PRINT "WRITING TO DISK --> "Fis (K) 

i S090 OP BM, "O".1,FIS(K) +": "+DNe 

1351906 Soelhierd i es EGR i=?) TG N +: PRINT #1,X¢I,K)3 : 
NEXT 





ian 10) CRG Seee- irae GOSUD S500 s IF EX=1 THEN 
Poi 30 

tat 20 FO=FO+1 =: FS(FO)=F1S(K) = GOSUB 3400 

Bad SO NEXT 

15140 RETURN 

199000 “ERROR HANDLING ROUTINE 

19015 IF ERR=10 AND ERL=8520 THEN NS="" : RESUME 8520 
feos) IF ERR <> 106 THEN. 19100 

freee ITF ERL=35100 THEN RESUME 3110 

19022 IF ERL=1040 THEN GOSUB 1100 : RESUME 1950 

ioe s IF ERL <> 12350 THEN 19100 

19030 IF SL=1715 COR SL=S010 OR SL=S4190 THEN PRINT “NO 
meee oeoN THIS DRIVE. FRESS ENTER 10 CONTINUE’: INPUT AS : 
ea-o =: RESUME 1300 

19043 IF SLIisAB240 OR SL=8330 OR SL=8550 THEN GOSUB 2200 : 
festiMe 1210 

Pee PRINT “"ERR="sERR; “ERL="sSERL =: ON ERROR GOTO 9 : 
peSsUuMe 


24 


ee ee pei eee ee eee 


10 CLEAR 1000 

20 ON ERROR GOTO 19000 

SO DEFSTR Z 

40 DIM FS#(935) .M$(11),N8 611) 

30 LAS=CHRS (123) : RES=CHRS (125) : MBS=CHRS (176) : BLSE=-CHRS (252) 
60 ZS="HFRHHS": ZZ="HH. BACCO" 

eo CLS 

BO FRINT TAB(20);"DATA FILE MANAGEMENT PROGRAM" 

90 PRINT:PRINT: PRINT "MENU: " 

100 PRINT "(€1) Create New Data Files” 

110 PRINT "(2) Edit Existing Files" 

mee FRENT *(S) Purge Existing Files" 

moo PRINT "(4) Print Data Files" 

140 PRINT "(3) FReturn to Master Menu" 

mae P8963 5=5:GO5SUB 190 

140 ON IO GOSUB 1000,1100, 1200, 1300.180 

7. SeT@ 7G 

180 RUN "“"MENU/BAS" 

190 REM k X® KEYBOARD INPUT SUBROUTINE £ 

200 PRINT oKPt+S, "OPTION DESIRED ";3LBS;MBS: RRS: 

210 ZI=INKEYS: IF ZI=""" THEN 210 ELSE IO=VAL(‘(ZTI) 

Meo PRINT DEP+21,ZI::FOR [=1 TO 100: NEXT 

Seem i821 OF T6255 THEN PRINT OEP, BLE3:PRINT OKP, "4 
NOT A VALID OPTION #*€84"53: FOR I=1 TO 1000: NEXT: FRINT 
wEP, BLS:;:GOTO 2ca0 

240 RETURN 

230 DNS=INKEYS: IF ONS="" THEN 250 ELSE T=VAL (DNS) 

250 IF [£0 OR I23 THEN 220 ELSE RETURN 

1000 REM kK * * NEW FILE ROUTINE * 

1010 CLS : INPUT "How many rows";RT 

1020 INFUT "How many columns";CT 

19050 BT=1 

1040 GOSUB 1409000 : GOSUB 1209000 : N=1 : CR=1 : GOSUB 46040 
Maa CLS:GOSUB 15050: RUN 

mee SOSUB 1500 : CR=1 : GOSUB 6040 

11190 GOSUB 12000 : RUN 

1200 GOSUB 3000 : RUN 

weed GOSUB 1500 : GOSUB 4000 : RUN 

1500 REM *® * * READ DATA ARRAY % * 

foto §86Diés="/" 

meee CLS PRINT “Disk Drive Containing Files (Q-3)%7" 
foeo GOSUB 250 

1540 SL=1715:GOSUB 2400: TF SL=0 THEN GOTO 13520 

maoo SOSUB 2300 : PRINT “" 

1360 PRINT: IF I10 = 2 PRINT “Enter the File To Be Edited": 
ELSE PRINT "Enter the File Tos Be Frinted"; 





p27 0 
1580 


INPUT NY 


IF NVSFO THEN PRINT:PRINT" £ &* * INVALID FILE # & 


aeesGOTO 1560 


1390 
1400 
1510 
1620 
1439 
1540 
1650 
1560 
14570 
1680 
1691) 
1790 
1710 
1720 
1730 
1740 
2OQ0 
2O16 
2020 
2OSO 
2O8D 
=OS0 
=060 
2070 
2OB0 
2090 
2100 
2110 
2400 
~410 
2420 
2430 
2440 
2430 
2460 
=470 
2480 
=490 
ZOO 
eal 
=B809 
=8190 
=820 
2330 
~B4o 
=850 


2850 


PRINT 


Fis=F$ (NV) 
OPEN "I",1,F1% 
INPUT #1.RT 
INFUT #1,CT 
INPUT #1,BT 
CLOSE 1 
GOSUB 12000 : 
OPEN "I" ,1,F1i% 
INPUT #1,RT: INPUT #1,CT:INFUT #1,BT 
FOR I=1 TO RT 
FOR J=1 TO CT 
INPUT #1,X¢I,J) 
NEXT J 
NEXT I 
CLOSE 1 
RETURN 
REM x * x CHECK 


PRINT 


Pemertee DPIRECTCRY 2 & 


IF FO=0 THEN RETURN 
P2=0 
FOR Pi=1i TQ FO 
Beis ime... S(Pli-Peo 
CEeQSse. i 
MEAT Pl 
Te SO iach [BRN ELSE GOSUEB 2400 : RETURN 
aoe 10 CORRECT DIRECTORY 
Wei t e—PO THEM 2110 
eae Seat es ESC Te HrFSCI+1) =: NEXT I 
be—PO-lee PS=P2tl =: RETURN 


REM * * * READ DIRECTORY OF FILES « + 
CLS 
F$="CTABLE/DIR: "+DNS 
OFEN "I",1,F% 
INPUT #1 ,F0 
IF FO=0 THEN 2490 
FOR I=1 TO FO 
INPUT #1,F%(1) 
NEXT I 
CLOSE i : 
CLOSE 1 
RETURN 
REM « x x 
CLS 
IF FO=0 THEN 2910 
FOR I=1 TO FO STEP 4 
FOR J=1 TO 4 
IF I+J-13>FO THEN 2910 
IF J=1 AND LEN(FS(I+J-1))<12 AND I>i 


GOSUB 2000 


Piorlay DIRECTORY @©€3% = 


THEN 





23/70 


Se ae (tose —i vy) s CS STRS(Ii+3-1)5") 


"SES (I+J—-1);3 


=sgo 
=~B890O 
29700 
2910 
SOO 
3210 
e270 
ee 

3240 
S200 
3250 
S270 
S400 
S419 
S420 
34359 
S440 
5450 
2460 
3470 
S480 
35500 
S610 
3620 
24650 
2640 
5650 
4000 
4010 
4020 
4230 
4040 
4950 


NEXT J 
NEXT I 
PRINT 
RETURN 
REM * * * CREATE 
DNS= an i) a9 
CLOSE 1 
F$="CTABLE/DIR: "+DNS 
OPEN "O",1,F% 
PRINT #1.0 
CLOSE 1 
RETURN 
REM «* x * WRITE NEW FILES IN DIRECTORY * * 
F$="CTABLE/DIR: "+DNS 
OPEN "O".1,F% 
PRINT #1,F0 
FOR I=1 TO Fa 
PRINT #1,F$(1)3","3 
NEXT I 
CLOSE 1 
RETURN 
REM * * * CHECK 
EX=0 
FOR I=1 TO FO 
IF FS=F$(1) THEN EX=1 : 
NEXT I 
RETURN 
REM * * &* PRINT ROUTINE x« x 
CLS : INPUT "Input Printer Width";FW 
PF=O0 : ML=-1 : DL$="'" : TL=-1 
FOR I=1 TO N 
FOR J=1 TO CT 
IF LEN(STRS(X(1,3))) 


mero e/ DIR” es xX 


IF FILENAME WAS ALREADY USED * X 


RETURN 


Pa ee 


ML=LEN(STRS (4 (1,9))) 


49060 
4070 
4080 
40990 


NEXT J 
NEXT I 
MC=Mi+1 
CLS: PRINT: PRINT 


"Enter Column Headers. Separating 


lines with" 


4100 


PRINT"the delimiter *'* to produce a two column 


header" 


41190 
4120 
4130 
4140 
4150 
4160 
4170 
4139 


FOR f=1 TO CT 
PRINT: PRINT “Enter the Column Header fer Column";!I 
INPUT MS(T) 
NEXT I 
MS=-1 
FOR Y=1 TO CT 
P1=0 
Pe tMoR (Pi+i, Mt (1). DLS) 


es 





41979 
4200 
4210 
4220 
42350 
4240 
4250 
4260 
4279 
4-280 
4290 
4500 
4316 
45209 
4330 


feet Dh (MS(1),Piti1,L) +STRINGS (SB, " 


4340 
A350 
4340 
4370 
480 
4390 
4400 
4410 
4470 
4430 


between columns is"3MS;". 


4440 


=P 2—-P 11 
IF P2=0 THEN L=LEN(MS$(1))-PF1 
IF L>MS THEN MS=L 
IF P2=0 THEN 4250 
Pi=P2 
GOTO 4180 
NEXT I 
IF ML>MS THEN MX=ML ELSE MX=MS 
FOR I=1 TO CT 
Pi=0Q : NL=0 : T$="" 
P2=INSTR(P1i+1,M$(1I).,DLS) 
L=P2—-Pi-1 
IF P2=0 THEN L=LEN(MS(1I))-Pi 
SF=INT ((MX-L)/2) : SB=MX-L-SF 
TS=TS+STRINGS (CSF, " 
a) 
NL=NL+1 = IF NL?TL THEN TL=NL 
IF P2=0 THEN MS(I)Y=TS =: GOTO 4380 
Pi=P2 
GOTO 4290 
NEXT I 


PRINT:PRINT "Would you like to mumber the rews (Y/N)? 
AS=INKEYS: IF AS=""" THEN 4400 ELSE PRINT AS 

he AS—"Y" THEN FPI=S ELSE PI=0 

MS=INTCCPW-MX x CT-PI)/CT) 

FRINT:FRINT"The maximum # of spaces that can be used 


PRINT: INFUT"Enter the number of spaces between 


columns”; SP 


4450 


N":TS 


4460 
4479 
4450 
4490 


te 2e-O THEN INPUT 
>: IF TS="Y" THEN 4020 ELSE RETURN 
fie oF. 0 OR SPSMS THEN RETURN 

Peers INMPFuUT “Enter Resort Title"; Ts 

eae 1 PRAM ATS) +PI+MX—-1 

mee (Pis2) =: P2=PW-P1 : 


PC hee telus T¥me Y ELSE 


TE LBNK TS) CPW THEN 


TS=LEFTS (TS, FW) 


4500 
4510 
4320 


"—"s oNEXT I: 


4530 
4540 
4550 


43560 


IF LEN(T$)3P2 THEN P1=INT ( (PW-LEN(TS$))/2) 
LPRINT TAB(P1)3;TS 
LPRINT TAB(PLi): : 
LPRINT 
LPRINT "" : LPRINT 
FOR J=1 TO TL 

FOR I=1 TO CT 

LPRINT 


SOR Si=1 16 ENTS) = LPRINGT 


fee T=) PAMKASP) +I) sMI DOS (CMS CI), CJ-1)X+1,MA4): 


43570 
43580 


NEAT I 
LPRINT 


4590 NEXT J 


4600 As="""— : 


Pome l—a TOAMK =< SS—ASt"—" =: WEXT I 


4610 FOR I=1 TO CT 





A520 EPRGAT TAB(CI-1) #(MXA+SP)4+P1) s Ass 

46350 NEXT I 

4540 LPRINT "“" os K=PI4+MX—-1—-INT ((MX—-ML) /2) 

4650 FOR I=1 TO RT 

4560 Peet —o THEN LPRINGT STRS(LIs "2 "3 

4679 FOR J=1 TO CT 

4680 LPRINT 

mes ((J-1) 24MX+SP) -LENC(STRS (X61. 3) ) 24K) XCILd) 5 
4590 NEXT J 

4700 DPR PSO tS 

4710 NEXT f 

4729 RETURN 

2000 REM * * * PURGE DATA FILE DIRECTORY & x 
gotO CLS 

9020 PRINT "Disk Drive Containing Files to Be Purced 
ae) 

eesO GOSUB 250 

70490 SL=5010 : GOSUB 2400 : IF SL=0 THEN GOTO 50190 
3050 K=0 

5960 PRINT TAB(20);" £EXkKX PURGE DATA FILES ¥xKx4" 
3970 FOR I=1 TO FO 

9080 AS=""SPRINT FSCI-K);" pote Lhe ey /N/E) "s 
sem. FZI=INKEYS:1F Zi="" THEN S090 

oem 1 7I=VON THEN FRINT 2i1:GOTG 3130 

Seeeeim 2i—"y¥° THEN PRINT ZYsGOSUB 5150 ELSE PRINT "N" 
3i2OQ NEXT I 

21iSO GOSUB 3400 

9140 RETURN 

geo Kitt &S(I—-K) 

moe tf I-Kk=FO THEN 35180 

3170 FOR J=I-K TO FO-1 : FSC(J)=FS¢(I4+1) +: NEXT J: 
3180 FO=FO-1 : K=K+1 

aeeO RETURN 

6000 REM &%xX EXAMINE DATA x 

S010 INPUT CR 

eee IF CR < 1 THEN CR=1 

5030 IF CR > N THEN CR=N 

ems0 ES=CHRS(S1) =: NS="" =F 1=0: CLS 

eea0 Cl—1:Ri=1:Bi=1 

6960 CR=R1:CC=C1 

6070 GOSUB 67650 

geel GOSUB 46240 : GOSUB 6299 


60990 Pi=0 

6100 AS=INKEYS =: IF AS="" THEN 6100 

6110 IF AS=CHRS (91) AND P1i=0 THEN NG="": GOSUB 5360 : GOTO 
6100 


6120 IF ASB=CHRS (27) AND P1=0 THEN Nt="":CR=1:GOSUB 
6380:GOTO 6100 

6130 IF AS=CHRS (10) AND P1=0 THEN N@="": GOSUB 6410 : GOTO 
61900 

6140 IF AS=CHRS (25) AND P1=0 THEN NS="": CR=RT: GOSUB 





6440: 


6150 
6100 
6140 
6100 
6170 


65005 


6180 


6560: 


6190 
6720 
6200 
6210 
6220 
ee 3 
6250 
6240 
5250 
6250 
6270 
6280 
6290 
6500 
6510 
6520 
E320 
5340 
63530 
6350 
6370 
5380 
63970 
6406 
6410 
6429 
6430 
5440 
64509 
6460 
6470 
6480 
6490 
6500 
6510 
6520 
65:30 
6540 
6550 
63560 
63570 
6380 


GOTO 4100 

IF AS=CHRS(8) AND P1=0 THEN NS="": GOSUB 5470 : GOTO 
IF AS=CHRS(9) AND P1=0 THEN N$="": GOSUB 6530 : GOTO 
IF AS=CHRS(24) AND P1=0 THEN N$="":CC=1:GOSUB 

GOTG 4100 

IF AS=CHRS(25) AND P1=0 THEN NS$="":CC=CT:GOSUB 

GOTO 4190 

IF (A$>"/" AND AS<":") GR AS="E" THEN Pi=1 : GOSUSB 
: GOTG 6100 

IF A$=CHRS(13) THEN Pi=0: GOSUB 4590 : GOTO 4100 

IF AS=CHR$(8) AND Pi=1 THEN GOSUB 6740 : GOTO 6100 
IF AS="D" AND Pi=0 THEN RETURN 

IF AS=CHR$(31) AND P1=0 THEN RUN 

GOTO 6100 

REM «x x x CLEAR ENTRY AREA «x x 

cc=1 

PRINT a2," " 

PRINT 3 2.CHRS(95) 

RETURN 

REM * * ® BRACKET ENTRY * x 


K=320+129% (CR-R1}+11*(CC-C1) +5 
TY=15360+K 
POKE TY,183 


POKE TY+11,187 


RETURN 

Reese tyY.,.t2o + PORE TY+i1l,128 
RETURN 

REM * * * PROCESS UP ARROW 4 &£ 


CR=CR-1:I1IF CReil THEN CR=1:GSCTO 6400 

IF CR<R1I THEN RIi=CR:RL=R1+4:6GO0SUB 57560 
GOSUB 6340:GOSUB 52990 

RETURN 

REM «x x xX 
CR=CR+1 

IF CRSRT THEN CR=RT:GOTG S460 

IF CRSRL THEN RL=CR:Ri=RL—-4:G0SUB 6340:GO0SUB 46750 
GOSUB 6340:GOSUR 46290 

RETURN 

FEM x xX # 
CC=CC-1 
teeweo.! THEN GE=1:GOTO 6520 


PROCESS DOWN ARROW # & 


PROEESS LEFT ARRQW mm & 


Beet TEN Cl=CEsCL=€17t4:GOSUB 6750 
GOSUB 6340:GOSUB 62909 

Re | URN 

Pen Gem @ PRECESS RIGHT ARROW <£)'4 
eC -Ce+ 1 

fea ee Ct TREN GOSUB 7000 


ie ee | SEMEL —E Cs C1=—CL—4:60SUB 6760 
GOSUB 62340:GOSUB 6290 


RETURN 


60 





6590 
5600 
6610 
66290 
6630 
5640 
6650 
6660 
5670 
5680 
5690 
6700 
6710 
ee 20 
6730 
CHRS ¢ 
5740 


6730 
6760 
CL=C1 
6770 
5780 
5790 
CHRS ¢ 
STRIN 
56800 
6810 
5829 
6830 
S840 
6850 
5860 
USING 
5879 
5880 
5890 
5900 
7QOQ0 
7019 
7020 
7030 
704G6 
7050 
7060 
11000 
11010 
11020 
11030 
11040 
11050 


REM * * & CHANGE VALUE IN FIELD x x 
GOSUB 6240 

IF N¢S="" THEN GOSUB 46530 : RETURN 

X (CR, CC) =VAL (NS) 

Né= as 08 

K=326+128% (CR-R1) +11% (CC-C1) 

IF X(CR,CC) >999999 THEN GOTO 6690 
PRINT 2K+1,CHRS (1298): 

PRINT USING Z$;X(CR,CC); 

RETURN 

K=K-2:PRINT OK,CHRS (128); 

PRINT USING ZZ;X(CR,CC); 

RETURN 

REM x x * PROCESS NUMBER ENTRY x # 
NS=N$+A$ : PRINT 92," "3s PRINT 9 2,N%;:F'RINT 
95) : RETURN 


IF LEN(NS) >O THEN NS=LEFT$(NS,LEN(NS)-1) : FRINT 2 2," 
": PRINT 2 2,N$ 

RETURN 

PRINT 2196, CHR$(143) : IF CT<C1+5 THEN CL=CT ELSE 


+4 

IF RT<R1i+4 THEN RL = RT ELSE RL=R1+4 

FOR I=Ci-1 TO CL-1 

PRINT 2197+11%(1+1-Ci) , STRINGS (4, CHRS (143) ) 3: PRINT 
128);:PRINT CHRS(4S5+1);:PRINT CHRS(128) ::FRINT 
G$(4,CHRS (143) ) 

NEXT I 

PRINT 2195,CHRS$(170):PRINT 2259, CHRS(170) 

FOR I=Ri-1 TO RL-1 

IF I > 8 THEN K1=0 ELSE Ki=1 

FRINT 2520+128k(1+1-R1)+K1,STRS(1+1);:PRINT CHRS(179) 
Sireed—Ci-—t We CL—1 

PRINT @527+11%(J3+1-C1)+128% (I1+1-F1) , CHRS (128) 3: PRINT 


7$:X (I+1,5+1) 

NEXT J 

PRINT 32387+128% (1+1-R1i).,CHRS$ (170) 

NEXT I 

RETURN 

REM * x * PROCESS LAST COLUMN ENTRY x 


Boe =: C1l=CE =: CL=C1+4 
CR=CR+1 <:- IF CR > RT THEN CR=RT 
IF CR > RL THEN RL=CR : R1=RL-4 
GOSUB 6240 
IF CT + S&S OR CR = RL THEN GOSUBR §760 
RETURN 
REM * * XK YERTIFY FILENSGME -—- FW=1 IF FAILS «& * 
FW=0 


P1=INSTR(1,5%,"/") 

IF P1 <> O THEN IF LEN(S$)-Pi23 THEN FW=1 : RETURN 
IF Pi >8 THEN FW=i : RETURN 

IF Pi=O THEN IF LEN(S$)>8 THEN FW=1 : RETURN 


ol 





11040 IF Pi 
ASC (MIDS (S$, P1i+1,1))>90 THEN FW=1 : 


<< > 


~~ 


O THEN IF ASC (MID$(S$,P1+1,1))<455 COR 
RETURN 


ie ir Sol (LENS (S5$.1)9<65 OCR ASC(LEFTS(S$,1))>90 THEN 


FW=1 RETURN 

moe) IF Pi=0 THENe P2=LENASS) ELSE P2=P1i-1 

11090 PS=1 

11100 FOR I=P3 TO P2 

111190 AS=MIDS(S$,1,1) 

ie i =O IF (ASC (AS) 747 AND ASC (AS)<358) GR (ASC (AS) +64 
AND ASC(AS)<91 ) THEN 11140 

mei 350 FW=1 : RETURN 

11140 NEXT I 

mera [F FPS=1 AND Fl <> O THEN PS=P1it+i +: P2=LEN(S$) =: GOTO 
111900 

11160 RETURN 

12000 REM kK * * DIMENSION ARRAYS x xX 

12010 NV=MEM-S00 : NV=INTC(NYV/ (44 (RT & CT & BT+1)))-1 

120206 IF MM > NV THEN PRINT "INSUFFICIENT SPACE TO READ IN. 
"sINPUT "PRESS ENTER TO CONTINUE"; AS 

meso DIM XC(RT,CT?) 

12040 MN=NV 

12050 RETURN 

15000 REM « K * RENAME FILES x x 

Peold CLS 

13020 PRINT "Disk Drive To Contain Files (0-3)%7" 

1350390 GOSUB 250 

135040 SL=8240:GOSUB 2400 : GOSUB 2800 : PRINT *"" 

1390350 INPUT “New Name";S$ =: GOSUB 11000 

15050 IF FW=1 THEN PRINT “INVALID FILENAME” : GOTO 
150350 

7 O F$=SS : GOSUB 3600 

13080 IF EX=1 THEN PRINT:PRINT "A File Currently Uses 
the New Filename. ":FRINT “Do You Want to Use it Anyway 
meny?” ELSE GOTO 135100 

15090 AS=INKEYS: IF AS="" THEN 15090 ELSE PRINT SS: IF AS< > 
mY" THEN 13050 

13100 Fis=Fs 

13110 RETURN 

14000 REM kK kK MAME FILE x x 

14010 CLS 

meos0 PRINT “Disk Drive To Contain File (0-3)7":GOSUB 250 
1405090 SL=83530 : GOSUB 2400 : GOSUB 2800 : PRINT "™” 

14040 PRINT 

14050 PRINT “FILENAME FOR TABLE"; : INPUT SS : GOSUB 11000 
14060 IF FW=1 THEN FRINT "INVALID FILENAME" : GOTO 14050 
14970 F$=SS$ : GOSUB 3600 

14080 Pee k— 1) Pee M PRINT YFILENANE IS ALREADY USED. TYPE 
Y TQ USE NAME ANYWAY ELSE N™s: INFUT AS : IF AS <> "Y" THEN 
14050 

140906 Fis=FS 

14100 RETURN 





aod 
150190 
eo 
T3050 
15040 
15050 
15060 
pe) 7 
15080 
aoOUSO 
Bat OO 
Poa 10 
Pol Zo 
io 1 SO 
151490 
moi SO 
melt 7 ©) 
LOO 
mat /O 
19000 


REM & & * EXITING ROUTINE FOR OPTIONS 1% 2% £ 


CLS 
PRINT "De You Wish To Rename File (Y/N)?" 
ZI=INKEYS:IF ZI="" GOTO 15030 
IF ZI="Y" THEN GOSUB 13000 
SL=8550 : GOSUB 2400 : GOSUB 2800 : PRINT "" 
PRINT: PRINT 

PRINT "WRITING TO DISK --> "SFiS 

OPEN "0".,1,F1$+":"+DNS 
PRINT #1,RT:PRINT #1,CT:PRINT#1, BT 


FOR I=1 TO RT 
FOR J=1 TO CT 
PRINT #1,X(1I,d); 
NEXT J 
NEXT I 

CLOSE 1 : 


SS—-Fis =: GOSUB S600 2 IF EX=1 THEN 


FO=FO+1 =: FS(FO)=F1S =: CGOSUB 3400 
RETURN 


REM x * * ERROR HANDLING ROUTINE x x 


19010 TF ERR=10 AND ERL=6590 THEN N$="" : RESUME 6630 
meee IF ERR <> 106 THEN 19080 

Meeso TF ERL=3150 THEN RESUME 3150 

Beeaqo IF ERL=2040 THEN GOSUB 2080 : RESUME 2050 

meso IF ERL <> 2430 THEN [9080 

Beveo 1 St=17i5 OR SL=5010 OR SL=5410 THEN PRINT “NO 
meeeo UM THIS DRIVE. PRESS ENTER TO CONTINUE”: INPUT AS 
SL=0 =: RESUME 2500 

19070 IF SL=8240 OR SL=8330 OR SL=85S530 THEN GOSUB 3200 
RESUME 2410 

meeeew PRINT “SRR="5;ERR3 “ERL="sERL =: ON ERROR GOTO © : 
fe oUME 





PmROGRArP 4. WITtLCOA SEAS 


Samm mmm mmm mE 


10 CLEAR 1000 

20 ON ERROR GOTO 19000 

Soe iN P:DEFSTR z:DEFDBL F,L 

femeoete SCLSO00, 1) ,F$(9S) ,P (1500) ,X(1500) ,BtS0) 

29 LBS=CHRS (125) : RBS=CHRS (125) : MBS=CHRS (176) : BLS=CHRS (252) 
50 REM * * * WILCGXSON SIGNED RANES TEST x &# 

79 CLS:PRINT 920, "WILCOXSON SIGNED RANKS TEST" 

72a PRINT:FPRINT: PRINT “Disk Drive Containing Data Files 
me)?" 

SO DNS=INKEYS:1F ONS="" THEN 80 ELSE I[=VAL (DNS) 

meee 120 OR [23 THEN 75 

Seeol—-100 =: GOSUB 11000 : IF St=0 THEN 70 ELSE GOSUB 12000 
seers iN; ** 

Meee RINT: INFUT “Enter the Two Data Files (e.g. 3,5) "3F1,F2 
fem FiceO OR F2>FO OR Fl<1 OR F241 THEN 100 

BeOOPEaN “I",1,FS(Fil)+"%: "+DNS 

150 INPUT #1.N 

140 FOR [I=1 TO N : INPUT #1,5°¢(1,0) : NEXT I 

fo CLOSE 1 

mee eoeemwe 1", 1 FS(F2) +": "+DNS 

i7O INPUT #1,™ 

Poor OoR I-11 FO M =: INPUT #1,S5¢4(1,1) : NEXT I 

Meo CLOSE 1 

200 IF N=™ GOTO 250 

2190 PRINT:PRINT "The Data Sets are NOT of Equal Sizes. Do 
You Want to Continue using the Smaller # of data points, 
Enter other Data Files. or Return to the Master Menu 
fee a7 RR)?" s 

Meo 2ISINEEYS:IF ZI=""* THEN 220 

feet 2 i="C” THEN PRIN? Z2i:GOTO 240:ELSE IF ZI="E" THEN 
meee eS UC SRUN ELSE IF Zi"R" THEN PRINT "R":RUN"MENU/BAS” 
eese GOTO 220 

=40 ITF N > M™ THEN N=M 


=20 J=1 

~50 FOR I=1 TON 

270 So. 1 )—S(inO)—-Sit,1) 
~B0 IF St(J,1)=0 THEN J=J-1 
pees ©) J=J+1 

500 NEXT I 

S10 N=J-1 


329 FOR I=1 TO N : X(I)=ABS(S(2T.1)) =: NEXT I 
SSO M=1 2: Ni=N Sone eoNives: Ftly=l =; NEXT £§ : GOEVE 


10000 

540 JS=1 

So0 FOR K=2 TG N 

350 R=ABS(X(P (JS) )/X (PCE) )) 
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370 IF R < .999 OR R > 1.001 THEN JE=K-1 : GOSUB 14000 
: JS=K 

S80 NEXT EK 

390 JE=N : GOSUB 14000 


400 TR=0: TN=0 

410 FOR I=1 TON 
420 fees Vl. ioe 
TN=TN+ABS (S(1I,1)) 
430 NEXT I 

440 IF TP > TN THEN TS=TN ELSE TS=TP 

450 CLS:PRINT 20, "WILCOXSON SIGNED RANKS TEST”: PRINT: PRINT 
4650 PRINT "HYPOTHESIS NULL 

ALTERNATIVE" 
4709 PRINT " 
ea y)" 

480 PRINT " 
may)” 

490 PRINT " 
Sexy)“ 

300 PRINT: PRINT 
Respectively” 


O THEN TP=TP+S(I,1) ELSE 


Se SS a ECX) =/= 


CZ) EYWGAy a= ty) EG - 


C5) Ste = ERY) Ss 


"Where X and Y represent Data Sets 1 & 2 


310 JP=640:KS5=2:GOSUB 8000 

320 IF ITOQ=1 V=TS ELSE IF [O0=2 V=TP ELSE IF IO0=35 V=TN 

320 CLS:PRINTS20, "WILCOXSON SIGNED RANKS TEST" :PRINT: PRINT 
340 PRINT “Test Statistic = "3V 

24S PRINT "Sample Size = ";N 

3990 GOSUB 2900 

269 ON IO GCSUB 600,800, 1000 

379 PRINT o99°60,"FPress SPACE RAR when ready to return to 


main menu"; 
=S0 IF INKEYS$= 


Pan ocoTra soo ELSE RUN “MENU/ EAS" 


e7O REM *e* *& NULL HYPOTHESIS #1 Kk &# 


meee FRING “P-Value = "s:FRINT USING 
620 PRINT: INFUT"Enter Desired Aipha 
630 IF AL<9O OR AL?1 THEN PRINT "#4 
BETWEEN O & 1":GOTO 4990 

5490 PRINT:PRINT “For Alsha ="sAL;". 
eee iF AL«2 * FOX THEN PRINT 
eao PRINT " 
67) 
580 
B00 
310 
B20 
S30 IF ALO OR AL2+1 THEN PRINT 
BETWEEN O & 1"°:GOTO 620 

S40 PRINT: PRINT "For Alpha ="s3;AL;", 
850 IF AL FOX THEN PRINT "“accept"s 

S360 FRINT " 
870 FRINT 

S80 RETURN 


RETURN 
Sore taet SHULL BYPOTHESIS 
FRINT "P-Value = 


"KA 


55 


the hypothesis that EX) 
PRINT "-vse- the alternative that 


zt] 2 YX 
apie No sic 
PRINTSINPUT"Enter Desired Alpha 


the hypothesis that EX) 
Yeyo— the alternative that E(X) 


"H#H#H#S" 52 & FOX 
Pevel";4L 
FRROR *£4% ENTER A VLAUE 


3 
"accept"; ELSE FRINT “reject; 


ee 


SO af eet 


"HH #HHH" S FOX 
level "sAL 


FRROR ¥#2£%* ENTER A VLAUE 


Sloe eat "reject"; 


= E(y)" 
Ey) a2 





1000 REM * * xk NULL HYFOTHESIS #3 # & 

1910 PRINT "P-Value = "3:PRINT USING "#8. #3##°:FOX 
1020 PRINT: INPUT”"Enter Desired Aipha Level';AL 

1030 IF AL<O OR ALS1 THEN PRINT "“*** ERROR £*% ENTER G 
VLAUE BETWEEN O & 1":GOTO <2O 

1040 PRINT:PRINT "For Alpha ="3sAbLs", "S$ 

1050 IF a. <FOX THEN PRINT "accept"; ELSE PRINT "reject"; 
1960 PRINT " the hypothesis that E(X) => ECY)" 

1070 PRINT "-vs- the alternative that EX) < ECY)." 
1680 RETURN 

2000 REM * * & WILCOXSON DISTRIBUTION & x 

2020 PX=V: IN=N:FLAG=0 

2O30 IF VEIN CINt+t1)/2 OR Veo CR IN<1 THEN PRINT " * &* ERROR 
¥ #£":GOTO 100 

2040 IF VEIN CIN+13/4 THEN V=IN& CIN+1) /2-V-1:FLAG=1 
20350 IF V327 OR CINSS AND VIN) THEN GOSUB 2600 ELSE GCSUB 
2070 

2060 RETURN 

2070 IF INT(V)=V GOSUB 2110 ELSE GOTO 20965 

-O380 RETURN 

2090 V=EINT (CV) 2:GOSUB 2110: TFOX=FO%: V=V+1:GCGOSUB 

2110: FOX=(FOX+TFOX) /2: G0OX=1-FOX 

2-100 RETURN 

“110 11¢4¢2)=6:11¢63)=02: 11 (4)=0:71(5)=0:1146)=0 

={120 C=1 

“130 N=INT(.5% (-14+S0QR(61+8xV))) 

2140 FOR I=1i TO WN 

2150 Ti¢i)=0 

=160 1T11¢41)=11¢41)+1 

ery IF I<2 GOTO 2320 

21820 Ti¢é2S)=11¢1) 

eesO Li (2)=11(2)+1 

2200 IF ITe#3 GOTO 2320 

2210 11¢S5)=11¢(2) 

meoo [1(S)=11(3)+1 

Beso IF <4 GOTO 232 

2240 T1if4)=11 65) 

2eo0 T1¢4)=11(4)+4+1 

Seow IF [<5 GOTO 2320 

2e/O T1iS)=11 64) 

Peo 11(5)=11(35)+1 

e270 IF [<& GOTO 2320 

20900 11¢46)=11 5) 

25310 Tit6:=11¢(6)+1 

2520 S=11¢1)411¢02)4+11¢35) +11 64) 4+7211(05)4+1106) 

2c40 IF SV THEN GOTO 23690 

Peeee €—C+1:0N 2 GOTO 23550, 235490, 25350, 2510, 2519, 2500 
2s60 IF I>1i GOTO 2380 

23/0 T1i¢1)=IN:GOTQ 2550 

20900 FOR IL=1 TO I[I-1 

20770 S=0:J=I-IL 
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2400 
~4190 
2420 
24350 
2440 
2450 
2450 
24790 
2480 
=490 
2300 
mio 
ZeLo 
eet) 
~o49 
20 
2360 
Bay OD 
Be fS 
aueO 
2a 
2500 
2410 
2S20 
Zoo 
2O4o 
ZOHa0 
2665 
2479 
=630 
2h9O 
2700 


IF J-1<=0 THEN GOTO 24490 


FOR K=1 TO J-1 
S=S+I1i (kK) 
NEXT kK 


S=S+CIL+1)1¢(J)+¢CIL+1) 4 ¢IL+2)/2 


Le S-2V THEN 
bi Ci ti—ibLjy=in 


ON J GOTO 2540, 2530 


NEXT Ib 
TIC =INI GGTO 2550 


i 
i 
IF 
oe 
ee 
ia 


I1¢S5)<1IN 
Pita) < LS 
Ti ‘4)<1IN 
Pe) 1M 
i272 t=)< I 
Piel >< TN 


GOTO 


GOTO 
GOTO 
GOTO 
GOTO 
GOTO 
Seta 


2430 


2310 
2280 
2250 
ee) 
2190 
~1&60 


95> 


te ee ee 


NEXT I 

FOX=C/2C IN: QOX=1-FOX 

IF FLAG=1 THEN FOX=00%:0@0X%=1-FOx%X 
RETURN 

REM * * * WILCOXSON SIGNED RANE 
MU=INS CIN+1) /4 

VYAR= (2N4+17) 4 CIN+LIN/ 24 

NX=V4+0.5 

GOSUB 7OQ190 

1 4=(SNC24+3N—-1)/7 C(1ON# CIN+1) % (2N41)) 

W= (NX—-MU) /SOR (VAR) sF=(EXP(-O0.S54WO2) ) /SOR(S. 2983198) 
PS=F Xk ONC S-S3RW) 

FOX=FOX+L4ASFS 

FS=—-FR(WLS-1LOAWCDCS+15%W) 

FYS=—-FR OWL 7-21 AWOS4+1OSFWOS-LOSsW) 


APPROX DISTIBUTION 


L6=4% (2N0C44+6N035—-SN41) / (S358 CINK CIN+1) XK CSN41)) 02) 


2710 
2720 
2730 
2740 
7000 
7910 
7020 
70zZ0 
7040 
7050 


FOX=FOXtL4é4F5+.5%F74L4C2 
IP FLAG=1 THEN FOX=1-FOX 


@OX=1-FOxX 

RETURN 

REM * * * NORMAL DISTRIBUTION #& «* 
SD=SOR (VAR) 


XN= (NX-MU)D /SD: AX=AHNS (XN) 
XZ=. S9R942S KEXP ((-AXSAX)D /2) 
NT=1/ (1+. 23164194AX) 


Gene 26 (1. SI SSSISANT)—(. SS65S5S58ENTC2)+(1.78147S94NT0S)-(1.921 
SCOANTE4)4+¢1. 550274 4NTOS) ) 


#060 
7979 
7080 
Pos) 
too 
7710 


IF XN<0Q GOTO 7990 
FOX=1-Q0xX 

RETURN 

FOX=Q0X 

BOX=1-FOX 

RETURN 
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8000 REM * * x 


CFPTION SELECTION ROUTINE & & 


BO1O PRINT JIP+tS. “NUMBER OF OPTION DESIRED 
mess CHRS (175) :RES 


BOZO ZI=INKEYS: 0F zr" 


THEM 8020 


BO030 PRINT d3P+31,2Z13: IQ=VAL(ZI) 


8040 IF TO<1i OR [O+KS 


THEN PRINT a2JP,"*%k NOT A VALID 


OPTIGQN &kx*"3:FOR T=1 TO 1000:NEXT:PRINT d09P,8LS:GOTO 8010 
10000 REM * * kK ARRAY SORT «x x 

10010 L=1 : BiL)=Ni+1 

190020 J=B(L) ’SET END OF ARRAY SEGMENT 

19030 T=M-1 “SET START GF ARRAY SEGMENT 

19040 IF J-M < 3 THEN 10236 "HANDLE 1 OR 2 ELEMENTS IN 
SFECIAL CASE 

10050 MI=INT(CI4+J) /2) Soe) GIP ARE ELEMENTS 

19060 REM FIND A LARGE ELEMENT AMONG THE SMALL CONES 
100796 =i+1i 

10980 IF JT=J THEN 101790 

109990 pee i Ce CdD). THEN 20070 

19109 REM FIND A SMALL ELEMENT AMONG THE LARGE ONES 
1O91id J=J-1 

mote) TF f=J THEN 10170 

19130 IF X(Pt(J)) S= XCPCM1)) THEN 10110 

19140 REM EXCHANGE OUT OF PLACE ELEMENTS 

19150 H=PCT) =: PCISHP(I) =: PCD)=H =: GOTO 190970 
101609 REM ADJUST COMPARE ELEMENT TO NEW MIDDLE 
SeeMeENT 

101790 fe 1 > Mi TREN i=i-1i 

10130 IF J=M1i THEN 10210 

190150 Bolo e eh) —P imi) -: PiM1)=H 

19200 REM SAVE STARTING POINT FOR ARRAY SEGMENT OF 
LARGE ELEMENTS 

moe 1 O L=L+i : BiL)=I 

moO2 SO GOTO 10020 

Meee IF J-M < 2 THEN 10270 

19240 Poe ee oe a — A CP PIt1) ) THEN 10260 

19250 Paes eto —e (litt) = P(M+1)=H 

19260 REM SET BEGIN AND ENDPOINTS FOR ARRAY SEGMENT GF 
miamee FL EMENTS 

mee7O M=Ht(L)+1 2: L=L-i 

10280 IF L>0 THEN GOTO 10020 

19290 RETURN 

fmeeoy REM £ * * READ FILE DIRECTORY * & 

119010 CLS: PRINTO2O, "WILCOXSGN SIGNED RANKS TEST": PRINT 
11020 FS="EDITDATA/DIF: "+DNs 

mes” GPEN "“I",1,FS% 

11040 INPUT #1,F90 

meee FOR f=1 TO FO 

11960 INPUT #1,°F$(¢1) 

11070 NEXT I 

meeo) CLOSE 1 

11090 RETURN 
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12000 
12610 
12020 
12030 
12040 
12050 
PRINT 
12060 


REM £ * * PRINT FILE DIRECTORY & & 
Sesto aw LEE GXSON SIGNED RANKS TEST": PRINT 
FOR I=1 TO FO STEP 4 
FOR J=i TO 4 
iP I+J-15FO THEN 12090 
IF J=1 AND LEN(FS(I+dJ-1))<12 AND T71 THEN 


PRINT TAB(16%(J-1))2"("2STRS(I4+J—-1) 2") 


"SFS(I+J-1); 


12070 
126089 
12090 


14000 FEM K *k & 


14010 


RETURN 


14020 


NEXT J 
NEXT I 
RETURN 
AVERAGE RANKS « & 
ee Oo JEneNM St twos et =SEN(S(PtUS) 13. ee US : 


AV=0 : FOR IT=JS TC JE =: AV=AV4+I = NEXT I : 


AV=AV / (JE-JS+1) 


14936 
sf 

14040 
19000 
19010 


Beeot—5010 OR SL=3S910) THEN INPUT 


Bae SS 
19100 


afm vomit ©. ot (1) 51) —-SEMASCP(1),1)3 4 AV : NEAT 
RETURN 

REM x * * ERROR HANDLING ROUTINE & XK 
IF ERR=106 AND {(SL=2005 OR SL=2410 OR 
Bie Sees ON hls 
RESUME 11080 
ON ERROR SOTO 0 


Se=o7 if GR Silo 
DRIVE. 
ENTER TO CONTINUE"SAS =: SL=0 : 


PRINT "ERR=":ERR:"ERL="sERL : 





— mc eS ee ee 


10 CLEAR 1000 

2O ON ERROR GOTO 199000 

Seer F:DEFSTR 2: DEFDBL FL 

40 DIM F$(95) ,P (S000) ,X (3000) ,B(50) 

30 LES=CHRS (1235) : RBES=CHAS (125) : MBS=CHRS (176) : BLS=CHAS (252) 
60 REM kK & XK MANN-WHITNEY TEST x X 

7O CLS: PRINT 220, "MANN-WHITNEY TEST" 

7o PRINT:PRINT:PRINT “Disk Drive Containing Data Files 

ro — 3 ae" 

SO DNS=INKEYS$: IF DNS=""" THEN 890 ELSE IT=VAL (DNS) 

meeiie iO Oe I>3 THEN GOTO 7S 

90 SL=100 : GOSUSB 11060 : IF SL=0 THEN 70 ELSE GOSUB 12000 
wee eeNT "" 

Boo ex PMTs IMPUT “Enter the Two Data Files (e.g. 3,53 "3F1,F2 
mo 0 Pi sFO GR F2SFO OR Fi<1 GR F2<1 THEN 100 

eee "l" terPStPid+"%: “+DNS 

130 INPUT #1,N 

140 FOR ~[=1 TON : INPUT #1,%¢(1I) : NEXT I 

1350 CLOSE 1 

aoe eee “I"*,1,FS(F2)+"3 "+DNS 

170 INPUT #1,™ 

180 FOR I[=N+1i TO N+M =: INPUT #1,.X¢2I) =: NEXT I 

aeo ClGSeE i 

200 MT=M 

210 Ni=N+M:M=1 : FOR T=1 TO Ni =: PtLI)=I : NEXT IT : GOSUB 
19000 

e20O JS=1 

250) FOR K=2 TO Nil 

240 TF XtP CIS) )<>X CP CK)) THEN JE=K-i: GOSUB 14000: UJS=k 
=J90 NEXT K 

260 JE=Ni : GOSUB 140900 

E72) 1S=0 

280 FOR I=1 TO WN 

mo 86 G=7TS+K (CT) 

SO0 NEXT I 

SOS UX=TS-NA(N+1)/2 : UY=NT-UX 

S10 IF UX < UY THEN UU=UX% ELSE UU=UY 

See CLS sPRINT 220, "MANN-WHITNEY TEST": PRINT: PRINT 


Seo PRINT “HYPOTHESIS NULL 

SiERNAaTIVE™ 

S40 PRINT " A) B(x) = EY? Etx)} =/= 
cay) ™ 

Sao PRINT * i SAA) eee) ECA? > 

a Yi" 

SS5OePRINT " (3) Etxy => (Esey) atx. (3 

eer) “ 





Ey fe) 
Resp 
380 
383 
THEN 
39) 
400 
405 
410 
429 
4350 
main 
440 
600 
610 
620 
630 
BETW 
540 
6350 
560 
679 
680 
800 
610 
B20 
830 
RE TW 
R40 
SS 
£250 
3790 
S80 
1000 
19010 
1020 
105 
VALU 
1040 
1050 
10S4 
1076 
1086 
2000 
2019 
2020 
2OSO 
2040 
20350 
ZO60 
xe 3G 


PRINT: PRINT "Where X and Y represent Data Sets 1 & 2 
ectively" 

JP=640:KS=3:CGOSUB 8000 

IF YO=i THEN U=UU ELSE IF I0Q=2 THEN U=UX ELSE IF I0=3 
U=uUyY 

CLS: PRINT920, "MANN: -WHITNEY TEST":PRINT:FRINT 

PRINT "Test Statistic = "5U 

PRINT "Sample Sizes are ";3Ns:PRINT " and "sMT 

GGSUB 2000 

ON I0 GOSUB 600, 800, 19000 

PRINT 29960, "Press SPACE BAR when ready to return to 
nenu"s 


IF INKEYS= "" GOTO 440 ELSE RUN "MENU/S5AS" 
REM * KX * NULL HYPOTHESIS #1 x xX 
PRINT "P-Value = "3:PRINT USING "#.#8##8#"52 £ FOX 


PRINT: INPUT"“Enter Desired Alasha Level ";4AL 

IF AL“<O OR AL?1 THEN PRINT "#4 ERROR £242 ENTER A VLAUE 
SeNeO & 1 scOTo! 620 

ERIN? sPRINT “Por Alpha ="sAL3", “$s 

Ir AL<S *% FOX THEN PRINT “accept”; ELSE FRINT "reject" 
FRINT "“" the hypothesis that E(X%) = EY)" 

PRINT "-vs- the alternative that E(x) =/= E¢(Y)." 
RETURN 

Bore Xe ee! MYPOTHESIS #2 «x 

PRINT "P-Value = "s:FRINT USING "“#. ####"5 FOX 

FRINT: INFUT"Enter Desired Alpha Level ":AL 

IF AL<O OR ALi THEN PRINT "4x ERROR £*X ENTER A VLAUE 
SEN O & 1°sGOTCG 620 

BaeNiseRIN? Sror Alpha ="sALs", "3 

Peeseoe ok tem PRINT “accept": ELSE FPRINT "reject"; 
PRINT " the hynothesis that ECX) <= ECY)" 


a6 


Print “-vs-— the alternative that E(x) > EtCY)." 
RETURN 

REM Oe es *© NULL HYPOTHESIS #3 x 8 

PRINT "P-Value = "3:PRINT USING "#.#R##8"3 504 


PRINT: INFUT"Enter Desired Alpha Levei";AL 

IF 4L<40 OR ALSi1 THEN PRINT "£24% ERROR £43 ENTER & 
Secerweeh Oo & 2*s6GOTG 620 

Pets PREMT Sor Alpha ="sALs"“, "3 

Poe He ae toe N FRINT “accept"™: ELSE PRINT "reject"; 
Brett. tine hypothesis that EX) => ECY)" 

PRINT "-vs- the alternative that E‘x? < E¢Y)." 
RETURN 

REM & * &* MANN-WHITNEY DISTRIBUTION x & 

M=MT 

IF M?N THEN M=N: N=MT 

FLAG=0 

tao Ny e—oS THEN U=N-U-1:FLAG=1 

IN=N+M: V=TS: [= 

IF USN GR U<O OR Net GOR Mei THEN PRINT " x x ERROR 
OTO i090 


Tae 





2079 
=~OBO 
2090 
=100 
BeALO 


2130: 


2120 
2130 
2140 
2150 
2-150 
2170 
2180 
220 
220 
2210 


2230 
2240 
2250 
2240 


2270 


ee F 
22.80 
a2 FO 


2500 


ere O 


S520 


279K) 


Sophie aareey 
pape Da | 


2o4td 
2o00 
23560 
70 
2Oo80 
2u9GO 
S400 
2410 
ZAZO 
2450 
T4A4ad 
S250 
2460 
2470 
~480 
7490 
2300 
2o10) 
ewe O 
2.) 
2540 
Coa 


23560 


ip Vee2e7 GR Meo THEN GOSUB 2500 ELSE GOSUB 2090 
RETURN 
teeter vy —Y THEM GESUB 2150 ELSE GOTO 2110 


RETURN 
V=INT(V):GOSUB 2130: TFOX=FOX:V=V+1:G0SUB 
FOX=(FOX+TFOX) /2:Q@0X=1-FOX 

RETURN 
11(2)=0:11 (3) =0:11(4)=0:11 (5) =0:11(6) =0 

C=0 

T1(1)=0 

Tici)=11(1)+1 

IF I<2 GOTO 2310 

11 (2)=I1 (1) 

11i(2)=11(2)+1 

IF I<3 GOTO 2320 

Psy =11 (2) 

Ti(S)=11(3) +1 

IF I1<4 GOTO 2320 

1164) =11(3) 

11(4)=11(4) +1 

IF I<5 GOTO 2320 

Ti (5) =11 (4) 

Ti(S)=11(5) +1 

IF 1<6 GOTO 2320 

11 (5) =11 65) 

11(S)=11 (4) +1 
S=I1(1)411(2)41163) 41164) 41105) 411 (5) 

IF S>V THEN GOTO 2350 

C=C+1:0N I GOTO 2540, 2530, 2520, 2510, 2500, 2490 
IF I>1 GOTO 2370 
Ii(1)=IN:GOTO 2540 
FOR IL=1 TO I-1 
S=0:J=I-IL 

IF J-1<=0 THEN GOTO 
FOR K=1 TO J-1 
S=S+11 (Kk) 

NEXT K 

eS eee Ce C1144) STL +2) /? 

IF S>V THEN GOTO 2470 

Ii (I+i1—-IL)=IN 

ON J GOTO 2530, 275270, 2510, 2500, 2490 
NEXT IL 
li¢i)=IN:GOTO 


2450 


=I40 


Bee li tS) < BM 
Pe 11(35)< Ee 
TF Lic4)<IN 
me F1t(S)< 18 
ce rel C27. DN 
eee t= TN 
S=1 


CECTOC 
GOTO 
GOTO 
GOTO 
GOTO 
GOTO 


2310 
2280 
2750 


2320 


le al ee 


ap 
2160 


POR r= TO N-1:S=St(IN-K) /(N-K):+NEAT & 


~1 
bd 





Zo? 
aot 
eae GU 
2600 
BOLO 
mow) 
mooO 
2540 


FOX=C/S: QOX=1-FOx 

IF FLAG=1 THEN FOX=00X: GOX=1—-FOXx 

RETURN 

REM K * kK MANN-WHITNEY APPROXIMATE DISTRIBUTION & xX 
MU=N4M/2: VAR=NEMS (N+M41) /12:3NX=U+.5 

GOSUB 7919 

L4= (MC2+NO2+IN+MAN) / C2OXMANS CIN+1) } 


LS=(2% (MO 44+N04) +4 4M4N48 (MO 24+N02) +64 (MOD) 4 OND 2S) +44 OME S4NO35) +7 
¥MANX COM+N) +MOC24+NO24+24M4N—-IN) / C21 0X (MODS) ONES) KC CIN+1) 02) 


2630 W={NX—-MU) /SGR (VAR) :F=(CEXP(-O. S¥WO2)) /SQR(S. 298319) 
2660 FS=F% (WOS-SHW) 

26790 FS=—-F# (WOS-LO*NCS+154W) 

eee F7=—-FPSé (WO7-21 AWE S+1OS4W0S-10S5kW) 
mo?! FOX=FOXt+L4 ¥€ FS 

a7 VO FOXS=FPGAtLS ¥£°FS 

mero) FOXA=FOX+.3 * F7 * L4C2 

2720 IF FLAG=1 THEN FOX=1-FQX 

2/59 Q@OX=1-FOX 

2/40 RETURN 

7000 REM * * * NORMAL DISTRIBUTION * & 
7O10 SD=SGR (VAR) 

7O20 XN=(ONX—-MU) /SD: AX=ABSCXN) 

7OSO XZ=. SI9BS42SKEXP (C C-AX#RAX) /2) 

7O4O NT=1/ Cit. 2PSL1E419KAX) 

7959 


Stee ee SIF SE ISMN) —C. SS6SO58S*NT C2) +01. 7814784NTC2)-(1.821 
ZOOANTO4)4¢01. SS502744NTCS) } 


meee IF XNO GOTO 7090 

mee FOxX=1-20X 

7080 RETURN 

49090 FOX=CO0X% 

meee CLOX=1-FOX 

71190 RETURN 

See eRe ketekx OPTION SELECTION ROUTINE * & 

Pere PRINT DJIP+S, "NUMBER OCF OPTION DESIRED 
ees CHRS (175) sRS8S 

S050 ZI=INKEYS:IF ZI="" THEN S020 

ee PRINT gJP+31.Z13: IO=VAL(ZI) 

S040 IF [0<1 GR IO+KS THEN PRINT oO0P,"xX*£* NOT A VALID 
BreeerOoM  *824"53FOR I=1 TO 1OG0:NEXT:PRINT dJP,BL$:GOTC 8010 
10000 REM & * k ARRAY SCRT «& x 

meotO tL=1 = BiL)=Nitl 

10620 J=B(L) *“SET END OF ARRAY SEGMENT 

100350 I=M-1 “SET START OF ARRAY SEGMENT 

10040 IF J-M < 3 THEN 10230 "HANDLE 1 GR 2 ELEMENTS IN 


meeelet CASE 


10050 MI=INT (CI +0) /2) (o2 COMPARE ELEMENTS 
19060 REM FIND & LARGE ELEMENT AMONG THE SMALL ONES 
10079 T=I+1i 

10080 IF f=J THEN 10170 


wy J 
(J 





10090 
10100 
190110 
19120 
10130 
10140 
101350 
19160 


eee 610) <= A(R CM1)) THEN, 10070 
REM FIND A SMALL ELEMENT AMONG THE LARGE ONES 
J=J-1 
IF {=J THEN 10170 
Peet) ee — ACP CMI) ) THERP 10110 
REM EXCHANGE OUT OF PLACE ELEMENTS 
Beet) wer C1) =F 0d t PCJO=H : GETO 10076 
REM ADJUST COMPARE ELEMENT TO NEW MIDDLE 


SEEMENT 


19170 
10180 
19190 
109200 
LARGE 
10216 
10220 
LOLS SO 
192740 
10250 
10250 
1 ARGE 
10270 
10280 
10296 
11000 
119190 
11620 
11030 
11540 


11050 


a 


ee ee ee 


ee) O 
22030 
12040 
12030 
PRINT 
12060 


IF I > Mi THEN I=I-1 
IF J=M1 THEN 102710 
H=P(1I) +: P{I)=P(M1) +: P(M1)=H 
REM SAVE STARTING POINT FOR ARRAY SEGMENT CF 
ELEMENTS 
L=L+1 : BtL)=! 
GOTO 10020 
IF J-M < 2 THEN 10270 
IF X(P(M)) X (P(M4+1)) THEN 10260 
H=P(M) : P(M)=P(M41) = P(M4+1)=H 
REM SET BEGIN AND ENDPOINTS FOR ARRAY 
ELEMENTS 
M=B(L)+1 : L=L-1 
IF L>O THEN GOTG 10020 
RETURN 
REM x * * READ FILE DIRECTORY * x 
CLS:PRINTS20, "MANN-WHITNEY TEST":PRINT 
F$="EDITDATA/DIR: "+DNS 
OPEN "I",1,F% 
INPUT #1,F0 
FOR I=1 TO FO 
INPUT #1,F$(1) 
NEXT I 
CLOSE 1 
RETURN 
REM * * * PRINT FILE DIRECTORY * x 
CLS: PRINT92O, "MANN-WHITNEY TEST": RINT 
FOR I=1 TO FO STEP 4 
FOR J=1 To 4 
IF 1+3-1>FO THEN 12090 
IF J=1 AND LEN(FS(I+3-1))<12 AND 


<= 
~ 


SEGMENT OF 


tee, ES 


Peemenes GlLomie— lt) ye Cs SIRS (1l4+g-1)3%") 


"SFSCI+J-1); 


12070 
12080 
129090 
14000 
14019 
14929 
149:30 
14049 


NMEXT J 
NEXT I 
RETURN 
REM * * X AVERAGE RANKS * £ 
iF JS=JE THEN X(P(JS))=JS =: RETURN 
AV=O:FOR 1=JS TO JE: AV=AVA+ISNEAT 1: 
FOR I=JS TO JE: X(P (I) )=AV: NEXT 1 
RETURN 


AV=AV/ (JE-JS+1) 


74 





19000 “ERROR HANDLING 

19010 IF ERR=106 AND (SL=2005 OR SL=2410 OR SL=2710 OR 3110 
OR SL=3510 OR SL=2910) THEN INPUT "NO FILES ON THIS DRIVE. 
PRESS ENTER TO CONTINUE"3AS =: SL=O0O : RESUME 11080 


i97100 PRINT "“ERR="S;ERR; "ERL="SERL =: ON ERROR GOTO 0 


“I 
Ci 





—FPROoOGRAr! ©. 


10 

20 

sO 

40 

30 

50 

79 

Lo—s 
eo 

oO 

160 
= iS 
110 
120 
Poe0 
140 
150 
1696 
179 
180 
190 
ZOO 
21 
22 
Eel) 
240 
242 
244 
2435 
246 
248 
“uO 
240 
ee () 
=2O 
290 
300 
J=J 
4 0) 


Sims 
wt 


380 
Es. 
S40 
co) 
SEO 
3/70 
So 


SPirtRnoOVv,’ BAS 


ON ERROR GOTO 19000 
Perot 2: DEFINT P 


eee 7 ye OO), Al (1500), 8150), X2(1500) ,x% (1500) 


Lt RS=CHRS (125) > RBS= CHRS (125) :MBS=CHRS (174) : BL$= BARS << 252 
REMee EO X SMIRNCOY TEST & & 
Ces he 20, MSOMIRNOV TEST” 


ES ONTsPRINTSPRINT 


ea 


DNS=INKEYS: IF 


“Disk Drive Containing Data Files 


DNS="" THEN 80 ELSE T=VAL (DNS) 


IF 40 QR I23 THEN GOTO 79 

St=1i0 : GOSUEB 11000 : IF SL=0 THEN 60 ELSE GOSUSB 13000 
RINT "" 

Siew se (NRW wemter the Two Data Files te.g. 3,5) "3F1,F2 
freee ischO 98 Fe-FO OR Fi<1 GR Fe<i THEN 110 

Beem "2 71 Fear i)t+":s °+DNS 


INPUT #1,N 


FOR I=1 TO N : INPUT 31,X¢1I) =: NEXT I 

CLOSE i 

NI=N-eM=i:sFOR I=1 TO Ne PCI =ItNEXT I1:GOSUB 10000 
FOR I=1 TO N:X1¢(I)=X(PC(I)):¢NEXT I 

meme: 2 oe Re Str Sl) +s “+ DNS 

INPUT #1,.™¢M 

Bee i=! TC MM =: INPUT #1.X¢1) : NEXT I 

Geese i 


MI=MMsMSITsFOGOR I=1 10 MMsPt(IJ=L2NEXT IsGCSUB 10000 
eee I—1 10 MMsAS (IV SHXCP CL) 2 sNEXT I 


IF MM="N THEN GOTO 2350 

mimeet—i TO NeX (I) =X1C1I) = NEXT I 

meme tal 10 MM s XItT)=H=X2CI) s NEXT I 

FOR IT=1 TO N =: ASCTI)=XtI) =: NEXT I 

T=N : N=MM : MM=I 

I=1:J=1: ACDOF=0: YCDF=0 

TP=0: TM=0 

IF IcN AND J?>MM THEN GOTG S80 

Poet ty AS CD) THEM I=It1isJ=Jti:sGOTO 33¢ 
Poteet AS CI) THEN GOTG SOO ELSE GOTO 3190 
Poem 1MSh I=I+isGOlTO SSO:ELSE IF Ji=MM THEN 
meee SSOsELSE GOTO Sao 

Pee ee) 6UTREN CETG S20 ELSE GOTO S30 
Sees iM THEN J=J+1 ELSE IF I <=N THEN JI=I+1 ELSE GOTO 


meuri=ti-1)/N 
p= ye DE—ACDF 


ACDF=(J—-1) /MM 


ee D< TM THEN 
fe DTP THEN 
SBan0 270 

SES TSRINTI2ZO, 


TM=D 
TP=D 


TSrIRNGy TES?’ SPRINT =sPRINT 





Seer eiNt “RYFGTHESIS NULL 
ALTERNATIVE" 

400 PRINT " ol} E(x) 
Gix)" 

410 PRINT " (2) P(X) = G(X) FOX) = 
S<x)” 

420 PRINT " (3) F(X) => GiX) Bix) 
ex)” 

430 PRINT:PRINT “Where FC(X) and G(X) represent Data Sets 1 
& 2 Respectively" 

440 JP=640:KS=3:GOSUB 8000 

pee  CeSsPRINTY2Z20, “SMIRNOV TEST":PRINT: PRINT 

450 ON I0 GOTO 470,480,499 

470 IF TP?ABS(TM) THEN D=TP ELSE D=ABS(TM):GOTO S00 

480 B=TP:GOTO Seo 

4990 D=ABS (TM) 

S00 PRINT “Test Statistic = “3D 

20S PRINT “Sample Sizes are "s3N3:PRINT " and "3M 

319 GOSUB 2000 

320 ON IG GOSUB 6900, 800, 1000 

oo0 PRINT 2960, "Press SPACE PAR when ready to return to 
main menu"s 


GX) F(X) =/= 


340 IF INKEY$S= "" GOTO 3540 ELSE RUN "“MENU/ BAS" 
S500 REM * *% * NULL HYPOTHESIS #1 8 & 
610 PRINT "P-Value = "3 :PRINT USING "B.8###8#"32 & 20% 


620 PRINT: INFUT"Enter Desired Alpha Level ";AL 

SO IF AL<G OR AL+1 THEN PRINT "xx ERROR «££ ENTER 4 VLAUE 
BETWEEN O & 1":GOTO 620 

640 FRINT: PRINT "For Alpha =" sAL:2". "Ss 

630 ITF AL<2400XK THEN FRINT "“accept"s ELSE PRINT "reject": 
6590 PRINT " the hypothesis that Fix) = Gtx)” 

670 PRINT "-vs- the alternative that F(X) =/= G(X)..." 

é80 RETURN 

200 REM * * * NULL HYPOTHESIS #2 «x x 

SiO PRINT "P-Value = "3:FPRINT USING "#. ####"; COX 

S20 PRINT: INFUT"Enter Desired Alpha Level";AL 

830 IF AL“9O OR AL?1 THEN PRINT "*3% ERROR #£4% ENTER & VLAUE 
BETWEEN © & 1":GOTO 46209 

840 PRINT:PRINT "For Alpha ="sALs", "3 

meee FL. OOX THEN PRINT “accept”; ELSE PRINT “reject”: 
860 PRINT " the hypothesis that F(X) <= Gtx)” 

BePePRINT “-vys-— the alternative that F(X) > G(X)." 

880 RETURN 

TOGO REM & & *& NULL HYPOTHESIS #3 « & 

1010 PRINT "P-Value = “s:PRINT USING "#. ####" 5 GOX 

1620 PRINT: INFUT"Enter Desired ASloha Level ";AL 

106350 IF AL<9O GR AL*S1 THEN PRINT "22 ERRGR **% ENTER A 
VLAUE BETWEEN © & 1":G60TO 620 

1040 FPRINT:PRINT “For Aipha ="sAL3:", "$s 

MOS IF AL<@0OX THEN PRINT “accept": ELSE PRINT "reject"s 
1050 PRINT " the hypothesis that F(X) => G(X)" 

1070 PRINT "-vs- the alternative that Ft(x) «< Gtx." 


ff 





1080 
2000 
2010 
2020 


FILSE IF (M/N)=INT (M/N) 


2030 

2040 

=O50 

2060 

2O79 

2080 

BOG0 

SO10 

oes 
BO20 

BOSO 

FO4o 

OPTIC 
10000 
1G010 
10020 
10030 
10040 
see l 
100306 
10060 
10070 
10080 
16090 
10100 
101190 
1O120 
190130 
10140 
109150 
19150 
eae e 
10170 
101980 
104199 
19200 
LARGE 
10210 
10220 
10230 
102406 
1023G 
10250 
LARGE 
10276 


RETURN 

REM *« * * SMIRNOV DISTRIBUTION xX & 

M=MM 

Pie eee oeme oti £2778 (. LP 2eSOR(N) +. 11/SQRiN) ) 


THEN S=D * S@R(N/ (N+M) ) 
oes ee SOR) ) ELSE S=D *teSQRt(N7 (NM) > 
faye & SUR (rT) ) 

GOX=(.135) 0(SC2) 

IF @OX<0 THEN GOX=0 

IF @OX?>1 THEN @0OxX=1 


FOX=1-00X 

PX=D 

RETURN 

Role see 2m CPTION SELECTIGN ROUTINE xX x 
Peo wewot tos WEBER OF OPTIGON DESIRED 


:CHRS (174): RES 
ZI=INKEYS:IF ZI="" THEN 8920 
PRINT 20P+31,21::I0=VAL(ZI) 
IF I0<1 OR IO>KS THEN PRINT 2J3P,"*x% NOT A VALID 
N  *«#k"3:FOR I=1 TO 1000:NEXT:PRINT 30F,BL$:GOTO 8019 
REM * * * ARRAY SORT x x 
L=1 : B(L)=Ni+1 
J=B(L) *SET END CF ARRAY SEGMENT 
I=M-1 *“SET START OF ARRAY SEGMENT 
IF J-M < 3 THEN 10230 ’HANDLE 1 OR 2 ELEMENTS IN 
AL CASE 
MI=INT ((1I49)/2) "SET COMPARE ELEMENTS 
REM FIND A LARGE ELEMENT AMONG THE SMALL ONES 
I=I+1 
IF I=J THEN 10170 
IF X(PC(I)) X(POM1L)) 
REM FIND A SMALL 
J=J-1 
IF I=J3 THEN 10170 
IF X(P(J)) X(P(M1)) THEN 10110 
REM EXCHANGE CUT OF PLACE ELEMENTS 
H=P(1) +: PCI) =P C9) P(J)=H : GOTO 19070 
REM ADJUST COMPARE ELEMENT TO NEW MIDDLE 


~ 


—_ 


¢ 
. — 


‘ 


TREN. 10070 


ELEMENT AMONG THE LARGE ONES 


- 


NT 
I - Mi THEN I=I-1 
J=M1 THEN 10210 
H=P CT) mee l—eMi) = PAPO SH 
Eeniesaveesot tae GING POINT FOR GRRAY SEGMENT Gr 
Se SeNis 
L=L+1 BiLj=I 
SHyea TOC2s 
ied «2 THEN 
IF X(P(M)) <= X(PCM4+1)) THEN 102860 
H=P (Mm) PAM =Pihrl)> P¢M+1)=H 
REM SET BEGIN AND ENDPOINTS FOR ARRAY SEGMENT 
SPEEMENTS 
M=BiL)+1i 


b-4 bd 


F 
F 


Toy 7 


OF 


=— 
— 


Bail 


oS 





10280 IF L>+O THEN GOTO 190020 

10290 RETURN 

11000 REM «x * * READ DATA FILE DIRECTORY x # 
11010 CLS:PRINTS20, "SMIRNOV TEST":PRINT 

11020 FS="EDITDATA/DIR: "+DNS 

11030 OPEN "I",1,F% 

11040 INPUT #41,F0 

11050 FOR I=1 TO FO 

11060 INFUT #1,F$(1) 

11070 NEXT I 

11080 CLOSE 1 

11090 RETURN 

13000 REM * * * PRINT DATA FILE DIRECTORY * «x 
13010 CLS: PRINTS20, "SMIRNOV TEST":PRINT 

13020 FOR I=1 TO FO STEP 4 


PeOSO FOR J=1 TO 4 

13040 Peet PO TimEM LS0FC 

ee SO TF J=1 AND LEN(CFS (14+ d-1))212 AND [+1 THEN 
aaoh?T 

123040 mee et mK CI 1 pee Ce SIRS C1tI—1) 3") 
eet itJ—i): 

13070 NEXT J 


13080 NEXT I 

13090 RETURN 

19000 *ERROR HANDLING 

19010 IF ERR=106 AND (SL=2005 OR SL=2410 OR SL=2710 GR 3110 
meee oon OR SE=S910) THEN INPUT "NO FILES ON THIS DRIVE. 
Beeso EMIER TO CONTINUE"SAS = SL=O0 : RESUME 11080 

Meee PRINT “ERR=“;sERR; “ERL="“sERL =: ON ERROR GOTO O 


Toe 





5S CLEAR 1000 

10 ON ERROR GOTO 19000 

moO DEFSTR Z: DEFINI P 

30 LBES=CHRS (123) : RBS=CHRS (125) : MBS=CHRS (176) : BLS=CHRS (S52) 
40 DIM F$(95) ,P(1500),X1(1500) ,B(50) ,X 1500) 

30 LES=CHRS (123) =: RBS=CHRS (125) : MBS=CHES (176) : BLS=CHRS (252) 
60 REM «x k x jLILLIEFOR TEST FOR NORMALITY « ¥ 

7O CLS:PRINT 220,"N ORM AL IT Y ee See 

SO PRINT:PRINT:SFRINT “Disk Drive Containing Data Files 

ae o)? "s 

90 DNS=INKEYS: IF DNS=""" THEN 90 ELSE FRINT DNS: I=VAL (DNS) 
100 IF I<0 GR I:35 THEN GOTO 80 

110 SL=110 : GOSUB 11000 : IF SL=0 THEN 70 ELSE GOSUB 12000 
werrmim, |" 

Pee PRINT: INPUT “Enter the Data File #"sFi 

130 IF FisFO OR Fi<i THEN 120 

Pee OreNM "I" 1, BS(F1id} +": "+DNS 

150 INPUT #1,N 

Mee FOR Ii=i TO N =: INPUT #1,4¢1) : NEXT I 

imo CLOSE 1 

mee CLSsPRINT 27458. “k *#€k ORDERING DATA 

5 x" 

190 Ni=N:M=1:FOR T=1 TO Ns PCIJ=HI2NEXT 1:GOSUR 10000 

200 SM=0: SUM=0: S2=0 

=<10 FOR I=1 TO 

ee 1 —X (P11) ) oe SMHSM+X1 (2) se OV=HX1 (1) -—SUM: SUM=SUM+BYV /1:S92=S2+ 
DVX¥(ZX1¢2I)-SUM) :NEXT I 

Seer elLSsPRINT "SELECT OFTION " 

Beo FRINT:PRINT:PRINT "¢1> MEAN and VARIANCE are bcth 
unknown" 

Seo FRINTSPRINT "¢2) MEAN is known, VARIANCE is unknown" 
eee PRINT:PRINT "¢3) MESN and VARIANCE are both known" 
=600 JP=896:KS=3:GCSUB 8On0 

290 ON 10 GOTO 300,320,540 

=00 MU=SM/N: VR=S2/ (N-1):SD=SOR(VR) - XCDF=9: T=0 

S1o GOTO 370 

329 CLS:INFUT “Enter the MEAN": MU 

250 VYVR=0:FOR I=1 TO N:VR=VR+¢X1¢1)—-MU) C2: NEXT 
T:SD=S@R(VRE/(N-1)):2:GGTOQ 370 

S49 CLS: INPUT "Enter the MEAN"; MU 

soo INPUT "Enter the VARIANCE"; VR: SD=SOR (VR) 

260 TF VR<O THEN PRINT "* x x* ERROR * * &* VARIANCE MUST 
eee OSITIVE":SOTO S50 

3/9 FOR I= 1 TON 

380 XN=(X1¢1)-MU) /SD: AX=ABS (XN) 

S70 XZ=. S9BS422XP 6 (-AXZAX) /23) 

4AOO NT=1/ 01+. 251641 98AX) 





410 


a0x= 


Zu x 
420 
430 
440 
450 
460 
470 
480 
4990 
SOO 
310 
3290 
St) 
San 
390 
S00 
BE TY 
610 
520 
630 
640 
BOD 
main 
310 
LOGO 
UNEN 
2020 
ZOO 
2040 
2OSO 
20560 
2O70 
=O0B9 
2090 
2200 
ZLOO 
2eOoO 
m2 1 


Be OD 
2210 
Zoog 
2510 
See, 
2530 
2U40 
2550 
2560 
ee 
2580 
2565 


eee) Fee om )— (6 GoGo oSCer NTE 2) +11. /B814784ENT03)-—¢C1. 822 


NTL4)+¢1.33527424NTE£S3 > 

IF XN<Q THEN @0xX=1-0Q0X 

FOX=1-00Xx 

D=AEBS (XCDF-FOX) 

IF D>T THEN T=D 

XCDF=I1/N 

D=ABS (XCDF-FOX) 

IF D>T THEN T=D 

NEXT If 

CLS:PRINTO20O,"N ORMAL IT Y 
PRENTSERINTSPRINT “Test Statistic 
PRINT "Sansle Size eh 

ON :9 GOSUB 2000, 2000, 4000 

FRINT “P—-Value "ssPRINT USING "#.88#3"; 004 
PRINT: INFUT“Enter Desired Alopha Level "3AL 
IF @AL<O OR AL?F1 THEN PRINT “44% ERROR X44 
Beheo oc 1"ss60To SIO 

PRINTS:FRINT “For Alpha ="3AL3s", "5 


eas a % 
al 


— 
— 


Beebe OOX THEN PRINT Yaccept"; ELSE PRING "reject"; 
PRINT " the hypothesis that the sample is from" 
PRINT “a normal pepulaticn."” 
PRINT 99609, "Press SPACE RAR when ready to return to 
menu's 

Poe re GEO SiG ELSE RUN "“MENU/BAS" 

REM kK X &* LILLIEFGR’S DISTRIBUTION — MEAN AND VARIANCE 
OWN 

DENCM=1 

D= T/DCENOM 

S=D% (—. 01+50R (N)+.285/SQR (N) ) 

TDENOM=CDENOM 

SESER 2200 

faeeon => .85 THEN DENGM = 1-—.24 * FOX:GOTO 2300 

emer. — 2s Sieh DENOM = i-Sa7 — ~.S6% FOX : GOTO 2200 
Poe — oo REN DENOGM = 1.3389 -— .é64% FOX : GOTO 
Pome — seo lise DENEM = 1.795 —- 1.12% FOX : SOTS 
BeNeM = 2.155 — 1.464% FOX 

TF ABS CTDENOM — DENOM) -991 THEN GOTO 26306 

RETURN 

Sexe ei ool (SL) 

TIF @OxX + 1 THEN ©0x = 1 

FOX = 1 —- OX 

RETURN 

Bieeee => -7O THEN. = ©. 74S4PR + 0.2799 : GOTO 2590 
IF PR => .30 THEN XQ = O.5464PR + 0.4187 : GOTO 2596 
Peer m => 620 THEM AG = 0.4644PR + 0.43588 =: GOTO 2290 
Peano ee time KO — O.48S4PR + 0.4525 : GOTO 2E96 
foots soe THEN XGo= 0.S5878PR + 0.4521 3: GOTC 2570 
20 = .S845KkPR + 9.4062 


Bi 


ENTER (p 


VLAUS 





eae t) 
2500 
25610 
S000 
VARIA 
20190 
SO20 
Or 
SOLS 
SO0SO 
SO4GQ 
2050 
SGSO 
SOTO 
oOo 
oe ld 
4000 
4010 
40920 
4030 
40460 
Rood 
B010 
aes 
S026 
BO 30 
8040 
ee nip 
BOSO 
10000 
19010 
10020 
1HOTO 
16040 


LEO9D 
19160 
19110 
LO 120 
10140 
19140 
19150 
190760 
ELEME 
10179 
1oidgod 
1o190 


10209 


XQ=XOQ/DENOM 
KQ=XO+.20/N 

RETURN 

PEM £st * LILETEFOR’S 
NCE UNENOWN : # & 
DENOM = 1: PX = 
DENOM PX /DENOM 
S=D* (.055+S0R (N) +. 293/SQR(N) ) 
TDENOM DENOM 

GOSUB 25090 

TET eOXx 
Ter OXx 


DISTRIBUTION —- MEAN ENOWN AND 


D 


- 
— 
== as 


= .9 THEN DENOM 
-95 THEN DENOM 
IF FOX fa9 oe DENGM 
DENOM pert ee 

IF ABS (CTDENOM—-CENCM) 
RETURN 

REM * * #£€ KOLMOGORCY DISTRIBUTICN * *# 
S=T* (0. 12+SQR(N)+.11/SQR(N)) 

OBOX=eIGol (SES) 

FOA=1-0% 


ioc — so tx 
1.47 -—- .4k FOX 
faa. sot Fox 


SEG S200 
SOTO 2200 
SORE S206 


~ 
- 

< 
~ 


3 
— 
< = 
~ 


by 


Olena GTO S020 


RETURN 
Bem © = & OPTION SELECTION ROUTINE & &% 
PRINT gJP+5, "NUMBER OF GPTION DESIRED 


SERRS (1745) RES 
ZI=INKEYS: IF ZI="" THEN 98020 
PRINT dJJP+S1,2Z13: L0=VALCZI) 
IF TO0¢<1 QR IO+KS THEN PRINT a3JUP,"*x* 
hi set's s POR 
RETURN 
REM 
at aes 
J=B(L) 
f=M— i 
Leo! 
CASE 
Mi=INTi¢CI4+d) /2) —Se?t COMPARE ELEMENTS 
REM FIND A LARGE ELEMENT AMONG THE SMALL 
I=I+i 
IF I=J THEN 
Boe Ct) 
REM FIND 
J=J-1 
T=J THEN 10170 
ie ce Cd) ee ely) ee EL 
Pe@eecenevGe OUT OF PLOECE ELEMENTS 
H=P iT) PCTI=P(J) = FCIJ=H GOTO 16970 
REM ADJUST COMPARE ELEMENT TO NEW MIDDLE 


NOT A VALID 
1=1578 {OCO0O!NEATIPRINYT gdUP.,. BLS: GOTO 3010 
kk * ARRAY SORT 
BR(L)=Nit+1 
“SET END OF ARRAY SEGMENT 
 Ss=eotan: CF ARRAY SEGMENT 
SoecreN 1220 “HANDLE 1 OR 


cod 


“~ 


x x 


2. Ra MEN awl 


SNES 


i} 
a” 


173 
poet (hid) THEN 100/70 
SMALL ELEMENT AMONG THE LARGE ONES 


1 
) 
a 


a 


< 


NT 
ie 1 
Gs 


Mi THEN I=I-1 

J=Mi THEN 102190 

H=P c(t) PCIJSP (M1) Si) 4 

REM SAVE STARTING FOINT FOR ARRAY 


a. 
al 


SEGMENT OF 





LARGE ELEMENTS 


meet © eae eee ia) = L 

10220 Gay 10020 

ceeco IF J-M < 2°tHEN 10270 

190240 tree teeeriye =— ACP KMart) ) THEN 10260 

10250 H=P(M) 2: PCM) =P CM+1) = & CM+1)=H 

10260 REM CSET BEGIN AND ENDPOINTS FOR ARRAY SEGMENT OF 


LARGE ELEMENTS 
10270 M=B(L)+1 : L=L-1 

10280 IF L>0O THEN GOTO 190020 

10290 RETURN 

11000 REM x * x READ DATA FILE DIRECTORY x x 

11010 CLS:PRINT920,"NORMALIT Y TE S T":PRINT 
110270 FS="EDITDATA/DIR: "+DNS 

11050 OPEN "I", 1.F3 

11040 INPUT #1.,F0 

11050 FOR I=1 TO FO 

11060 INFUT #1,F$(1) 

ive7O NEXT I 

11080 CLOSE 1 

11090 RETURN 

12000 REM «* x * LIST DATA FILE DI 
12010 CLS:PRINT9ZO,."N ORMALIT 
12020 FOR I=1 TO FO STEP 4 


PECTORY @ 'S 
ne Noe, ce ata eet Da 


12030 FOR J=1 To 4 
12040 IF I+J-12>FO THEN 12099 

12050 IF J=i AND LEN(F$(I+J3-1))<12 AND I>1i THEN 
PRINT =" 

12060 PRINT TAB(16%(J-1))3"(":STRS(CI+d—-1) 3") 

Ws FS(T4+I—1)3 

12070 NEXT J 


ees NEXT I 
£2090 RETURN 
19000 *“ERRCR HANDLING 


19010 IF ERR=106 AND (SbL=2005 OR SL=2410 OR SL=2710 GR 3110 
fee —e ot) GR SL=S910) THEN INPUT “NO FILES ON THIS DRIVE. 


fee = EMER TO CONTINUE™: 4S =: Si=0 : RESUME 11060 
meee PRINT “ERR="Z;ERR; “ERL="sERL : ON ERROR GOTO oO 





cme ee —— SR eee 


3 CLEAR 1000 

160 ON ERROR GOTO 19000 

EOmDErPSTR Z: DEFINT P 

30 LBS=CHRS (123) : RBS=CHRS (125) =: MBS=CHRS (176) : BLS=CHRS (252) 
409 DIM FS(9S) ,P (1500) ,X1(1500) , 8150) ,X (1509) 

30 LBES=CHRS (123) : RBES=CHRS (125) : MBS=CHRS (176) : BLS=CHRS (252) 
60 REM « * X* EXPONENTIAL DISTRIBUTION TEST * £ 

Peele Ss:PRINT 220,"E X PONENTIAL iE. Sain * 

SO PRINT: PRINT:PRINT "Disk Drive Containing Data Files 
aoa) COS 

0 DNS=INKEYS: IF DNS="" THEN 90 ELSE FRINT DNS: IT=VAL (DNS) 
100 TF [#0 QR I23 THEN GOTOQ 80O 

1190 SL=110 : GOSUB 11000 : IF SL=0 THEN 7O ELSE GOSUB 12000 
eereRiny “" 

120 PRINT: INPUT “Enter the Data File #";F1 

meet FizFO GR Fist THEN 120 

moemeeres TO 1s rPStriy+': *+DNS 

150 INPUT #1,N 

160 FOR IT=1 TON : INPUT #1,X¢1T) 2: NEXT I 

ime CLOSE 1 

BeOeeLlS:FRINT 2458, "* *€ kK CRODE RING DATA 

x x" 

190 Ni=N:M=1:FOR I=i1 TO Ne PCI) =H=I2NEXT 1:GOSUB 16000 

=00 SM=0 

210 FOR fT=1 TO N:exicilId=X(P CT) sSM=SM4+xX1i0TS sNEXT P:LA=N/SM 
Meee elS:PRINT “SELECT GPTION “* 

Meme PRINT:PRINT:PRINT “¢1) LAaMBDS is unknown" 

240 FRINTSFRINT "¢2} LAMBDA is known" 

2350 JP=896:KS=2:GOSUB 3000 

260 ON 10 GOTCG 290,270 

279 CLS: INPUT "Enter L&MEBDA":LA 

280 IF L& <= 0 THEN PRINT "* * ® ERROR *£ £ * LANEDA MUST 
Beer oolliIveE*:GOTO 270 

290 T=0: XCDF=0 

200 FOR I=1 TON 

Beer OX = L-EAPC-LALCI)) 

S20 D=ABS(XCDF-FOX) 

320 IF D > T THEN T=D 

=49 XCDF = I/N 

so0 D=ARS (XCDF-FOX) 

260 IF D > T THEN T=D 

eyo NEXT I 


Beemer S:PRINT@SO,."—E XP ONEN TIAL a Ee feo al Set 
270 PRINT: PRINT:PRINT "Test Statistic = "3T 

400 PRINT "Sampie Size = "3N 

4190 ON I0 GCGSUB 2609, 3000 

420 PRINT "F-Value = "s:PRINT USING ". ####": GOX 


84 





430 
440 
BE TW 
4350 
4&0 
470 
430 
490 
main 
BOO) 


2000 REM * * ¥* LILLIEFOR*S DISTRIBUTION — MEAN AND VARIANCE 


UNKEN 
2O2O 
2040 
ZO6D 
O80 
2100 
el2o 
2140 
2156 
~180 
22 OO 
fe -t) 


2250 
2280 
2460 
ZOO 
2460 
eo 
2469 
2O40O 
~460 
2060 
2469 
280 
=4&60 
2400 
24580 
2420 
2450 
~440 
2460 
2470 
SOOO 


S010 
5020 
SOSO 
5040 
SOS0 
5060 
2000 


PRINT: INPUT"Enter Desired Alpha Level ";AL 
ITF AL<O OR ALS1 THEN PRINT “k*x ERROR «kx 
FEN O & 1":GOTO 430 

Peden) Fer Alolba ="sALs;", “3s 

IF AL<QOX THEN PRINT "accent"; ELSE PRINT 
PRINT " the hypothesis that the sample is 
PRINT “a normal populaticn." 

PRINT 09960, "FPress SPACE BAR when ready to 
menu's 

IF INKEYS= 


"9? SOTO S00 ELSE RUN "MENU/BSAS" 


OWN 

DENOM=SOQR (N)+.25+.5/SQR CN) 

IF T * 1.3442/DENOM + .2/N THEN FOX = 1 : SOTO 2220 
FOX=.3:GOSUB 2250 

IF T}>XQ THEN Fi=FOX:F2=1:GOTO 2120 

Fi=0 : F2=FOQOX 

Ria = tr ithe) 7/2 

GOSUB 2280 

IF ABS(T-X¥O)<.0001 THEN GOTG 2220 

IF T>X@ THEN F1i=FOX ELSE F2=FQX 

GOTO 2120 

QOxX=1-FOQXx 

RETURN 

Te FPG@X => .95 THEN *0O = 3.125% FOA - 3.7808 =: GOTO 
IF FOX => .90 THEN XO = 2.026% FOX - 0.8357 : GOTC 
TF FOX => .80 THEN XQ = 1.124* FOX - 0.0249 : SOTO 
IF FOX => .70 THEN XQ = G.74Sx FOX + 0.2799 : GOTO 
TF FOX => .50 THEN XQ = 0.545%x FOX + 0.41937 : GSOTC 
IF FOX => .30 THEN xQ = 0.464% FOX + 09.4388 : GOTO 
IF FOX => .20 THEN XQ = 0.4858 FOX + 0.4525 : GOTO 
IF FOX => .10 THEN X@ = 0.587% FOX + 0.4321 : GOTO 
XO = .846*% FOX + 0.4062 

AQ=KOQ/DENOM + .2/N 

RETURN 

REM £ * kX LILLIEFOR’S DISTRIBUTION —- LAMBA KNOWN : 
S=T*(0.12+SGR(N)I+.11/SQR(N)) 

BOX=.1SEC (Ste) 

IF @OX>1 THEN ©0xX=1 

FOX=i1-Q0xX 

PX=T 

RETURN 

REM £ * *€ GFPTION SELECTION ROUTINE x x 


ENTER A VLAUE 


"reject"; 
from" 


return to 





8010 FRINT 2J3F+5, "NUMBER OF OPTION DESIRED 
"s1B$:CHRS (175) sRBS 

8070 ZI=INKEYS:IF ZI="" THEN 8020 

S030 FRINT 23P4+351,2Z1;:I1Q=VAL(ZI) 

8040 IF I0¢1 OR IO2KS THEN FRINT 30, "xxx 


NOT &A VALID 


CPt ON ekSS"5>2FGR I=1 FO 1O0O03NEXTS PRINT 2JP,8BL$:60TO 8010 
BoOoSO RETURN 

i90000 REM * * * ARRAY SCRT X* & 

10010 L=1 : BCLI=Ni+l1 

10020 J=Bi(L) *SET END OF ARRAY SEGMENT 

10030 I=M-1 “SET START OF ARRAY SEGMENT 

10046 IF J-M * 3 TREN 102350 “HANDLE 1 OR 2 ELEMENTS IN 


SeecitAl CASE 


10056 MiI=INT((I+9)/2) "SET COMPARE ELEMENTS 

10060 FEM FIND A LARGE ELEMENT AMONG THE SMALL ONES 
10070 I=I+1 

10080 IF I=J THEN 10170 

10096 IF X(P4I)) <= X¢POM1)) THEN 10070 

10100 REM FIND A SMALL ELEMENT AMONG THE LARGE ONES 
10110 j=J-1 

i9120 IF I=J THEN 10170 

10130 IF X€P6J)) >= XIP(M1)) THEN 10110 

10140 REM EXCHANGE OUT OF PLACE ELEMENTS 

10150 H=P (I) : PCID=P(I) +: PCJ)S=H : GOTO 10070 
19160 REM ADJUST COMPARE ELEMENT TO NEW MIDDLE 
ELEMENT 

10179 IF I > Mi THEN I=I-1 

10130 IF J=M1 THEN 10210 

10190 H=P (I) : PCI)SP (M1) +: PCMi)=H 

10200 REM SAVE STARTING FOINT FOR ARRAY SESMENT OF 
LARGE ELEMENTS 

10219 L=L+1 : Bi(L)=I 

10220 GOTO 10020 

10230 IF J-M < 2 THEN 10270 

10740 IF X¢POM)) <= X¢(P(M4+1)) THEN 10240 

10250 H=P(M) 2: POM) =P(M+1) : POM+1)=H 

19250 REM SET BEGIN AND ENDFOINTS FOR ARRAY SEGMENT OF 
LARGE ELEMENTS 

10270 M=B(L)+1 : L=L—-1 

10280 IF L>O THEN GOTO 10020 

10790 RETURN 

11000 REM +* * * READ DATA FILE DIRECTORY «* x 

11010 CLS:FPRINTS20,"—E X PONENTIAL TES T":PRINT 
11020 FS="EDITDATA/DIR: "+DN$ 

11050 OPEN "I",1,F% 

11040 INPUT #41,F9 

11050 FOR I=1 TO FO 

11060 INPUT #1,F%(T1) 

11070 NEXT I 

11080 CLOSE 1 

110906 RETURN 

12000 REM * * * LIST DATA FILE DIRECTORY x * 


846 





PeweeogemimbasO, “EE X PONENT IAL ee eomy smi 
feveoerOR I=1 TO FO STEP 4 


P2050 FOR J=1 TO 4 

12040 Reet ote THEN 1 2O9O 

re O50 IF J=1 AND LEN(FS(I+dJ-1))<12 AND I>1 THEN 
PRINT "" 

12040 pee et toe —1) ps MCS STRe(I+J-1)s:") 
mae (1 +J—1): 

12070 NEXT J 


12080 NEXT I 

P2070 RETURN 

19000 *ERRCR HANDLING 

Beer IF ERR=106 AND (SL=20035 OR SL=2410 GR SL=2710 OR 31190 
Sool —ooiO OR St=3910) THEN INPUT "NOG FILES ON THIS DRIVE. 
freoo ENTER TQ CONTINUE’: AS =: SL=0 : RESUME 11080 

meee o PRINT “ERR=";ERR; “ERL="sERL =: ON ERROR GOTO © 


CO 
“sd 





ery: em ee ee SSS 


3 CLEAR 1000 

19 ON ERROR GOTO 1909000 

mee DEFSTR Z: DEFINT P 

SO LBS=CHRS (123) : RES=CHRS (125) : MBS=CHRS (176) : BLE=CHRS (252) 
40 DIM FS(95) 

50 LBS=CHRS (123) : RBS=CHRS (125) =: MBS=CHRE (176) : SLS=CHRS (252) 
60 REM * * &* CHI-S@UARE TEST X x 

MeeeressPRINT @20,"C HI -SQG@UARE tT Eas 7" 

S80 PRINT: FRINT:FRINT "Disk Drive Containing Data Files 
os)? Ss 

909 DNS=INKEYS: IF DN¢="" THEN 90 ELSE FRINT DNS: I=VAL (DNS) 
meo FF 140 OR 23 THEN GOTO 80 

mee St=110 =: GOSUB 11000 : IF SL=0 THEN 70 ELSE GOSUB 12000 
meeeiNy 

£20 PRINT: INPUT "Enter the Data Fils #":F1 

Meo iF FistrFO OR Fi<1 THEN 120 

POrgreM “I".1,FS(F1i)+": "+DNs 

Boo INPUT #1,RT: INFUT #1,CT: INPUT #1,87T 

mee OihM X(RT,CT),RCRT) CCT) 

170 FOR Ri=1i TO RT 

mee FOR Ci=1 TO CT 

mo FOR Ri=i TC ST 

eee INPUT #1,.X¢(R1,C1) 

eee NEXT 51 

ge MEXT Ci 

moo NEXT Ri 

240 N=0O:0F=(RT-1)x* (CT-1) 

Gee FOR I=1 TO RT 

260 R(I)=0 

ee FOR J=i TO CT 

mew MCIS=R(T)+X (1,0) 

290 NEXT J 

S00 N=N+R(T) 

310 NEXT fi 

eee FOR J=1 TO CT 

S50 CfJI)=0 

540 FOR I=1 TO RT 

Beer (J)=CiJ)+X(1I,J) 

S60 MEXT I 
S70 NEXT J 
380 T=0:CZ= 
370 FOR [=1 
490 FOR J=1 TO CT 

410 E=R(1I)% ClJ)/N- 

Seem — 9 ((X(1,59)-E) £2) /7E 
eo e * SS THEN CF=CF+1 
440 IF © < 1 THEN CZ=CZ+1 


38 





450 
4&0 
465 
i. 
47Q 
4640 
490 
SOO 
er / ¢ 
mid) 
19 ; fe? 
ae 
ra ee 
S40 
SO 
BETW 
ae) 
a/ 
SAO 
=a9O 
500 
main 
510 
HOOD 
aS 
aHOGAd 
~OSO 
=~O8O 
aa A 31> 


air 


ee 


e200 
eee 
22 AQ 
ae 
ae 
Sit) 
gO 
SAO 
2200 
25g0O 
2406 
2420 
~44Q 
2460 
24EO 
SO0G 
— x 

SO1G 
5029 


NEAT J 
NEXT I 


hoa Nee h—< then T=(cAsS( (i, 1 eke X(2,2)-X 61,2) 
fee). can) hep * NY (RCL) SERS) CLL) Cl2)) 
Ses eivn a2, “ett | -~- SOUA RE wee)" 
PRINTS PRINT: PRINT “Test Statistic = "3sT 
FRINT "Degrees of Freedom ="; DF 
PRINT "Cells with expected frequency < S is "531900 
pees CTess PRINT "2" 
PRINT "Number cf cells with expected frequency < 1 is 
GOSUB 2000 
PRINT "P-Value = "s:FRINT USING ". ##8##"; 00% 
PRINT: INPUT"Enter Desired Aloha Level ";4L 

IF AL<G OR AL?1 THEN PRINT "#%4% ERROR *4£x ENTER A VLAUE 
FEN O & 1":GOTG S40 
BeemnsPRIMNt “Por Alpha ="sAL3", “s 

IF AL“<OCGX THEN PRINT "accent"; ELSE PRINT "reject"; 
PRINT “ the Khypothesis that the column cells” 
FRINT “have equal probabilities.” 
PRINT 2960,"Press SPACE BAR when ready to return to 
menu": 

Poeemieys— 9" COTO 610 ELSE RUN "MENU/BAS" 

Bee xk ® CHI-SQUARE DISTRISUTION «©. .% 

CA=T 

IF DF+l0o GOTO 2460 

DF=1 

eem I=INT(DF) 10 2 STEP -2 

DF=DP 

NEXT I 

ee erin 1 (DP +i) /2)) EXP (-CX/2) SDR. 

eee (DF /S)=DF/2 GOTO 2220s 

Geen C2 /CK/S.141593)* 4, 

GOTO 2240 ' 

}i=i 

{C=1 

MCc=1 

D=DF 

D=D+2 

MC=MC&CX/D 

TF MC<. 000001 GOTO 2400 

i C=LC+MC 

SOTO 2300 

FOX=Li4NP* LC 

Q@GX=i-F OX 

RETURN 

Peer Et iy Sy) — (1-(2/(9k 0F))}))/S0ORC2/ 09% DF) ) 
MU=05 : VAR=1 : [2=i1 

eM x * ¥ Mm O RPM A L Petes Gigix 2 2 uf I Gi 


SD=SER (VAR) 
AN= (NX-MU) /SD: AX=ABS (XN) 


a9 





S030 X7Z=.S989423% EXP ((-OXAK} /2) 

S040 NT=1/(1+.2516419%x AX) 

S050 

QOK=XZK(C. SISSBISENT) —¢€. SS6S56S3RNTOC2) +01. 7981478xXNTCIS)-(1.821 
SSOENTE4)4+¢1. 5S202744NTOS) ) 

S060 IF XN<0O GOTO 3090 

3G70 FOX=1-Q0X 

SO80 RETURN 

5090 FOX=00% 

21i0G@ OOX=i1-FOX 

5110 RETURN 

1272000 REM * * * READ DATA FILE DIRECTORY x & 
Meer CLS:FRINT 2290."C HI -S@G@uUARE TES T":PRINT 
11020 FS="CTABLE/DIR: "+DBNS 

most) CFEN “I",1,F% 

11040 INPUT #1,F9 

11050 FOR I=1 TO FO 

11050 INPUT #1,F3¢13 

11070 NEXT I 

11980 CLOSE 1 

11990 RETURN 


12000 REM «x * * PRINT OUT DIRECTORY CF DATA FILES «* x 
127010 CLS:PRINT 920,"C HI-SQ@QuUARE TES T":FRINT 
12020 FOR I=1 TG FO STEP 4 

12030 FOR J=1 TO 4 

12040 IF I+J3-1>FO THEN 12090 

17050 IF J=1 AND LEN(F$(1I4+J-1))<12 AND I?1 THEN 
PRINT 

12060 FRINT TAB(16*(3J-1))2"("3STRS(I+d-1)3") 
"SFS(I+d-1); 

12076 NEXT J 


12080 NEXT I 

Bee7) RETURN 

$9000 “ERROR HANDLING 

19010 IF ERR=106 AND (SL=2005 OR SL=24190 OR SL=2710 OF Fi 
SeeseaoolG GR SL=5910) THEN INPUT "NO FILES GN THIS DRIVE. 
memes Sli R TO CONTINUE’; AS = SL=0 : RESUME 11080 

Peers INT “ERR="sSERRs "ERL="sERL =: ON ERROR GOTO 0 





mee CREAR Soa0 

eeepeRinit I-K = DEFSTR Z 

SO eDEFDBL F,L 

40 LBS=CHRS (1235) = RBS=CHRS (125) =: MEBS=CHRS(175) : 
BLS=CHRS (2352 

30 CLS : PRINT TARB(2SO)s"**%e MENU 4K” 

SemereiNt™ (1) Frobabidity Distributions 

74 FPRINT' (2) Inverse Probability Distributions 

SO PRINT "(3) Return to Frocgram Selecticn Menu" 

GO KP=704 : ES=3 : GOSUB 130 

Po@eon 10 GOSUB 300, 1000,120 

moeeGOTO SO 

120 RUN "MENU/BAS" 

120 REM 

140 PRINT gdEP+S. "OPTION DESIRED ";LBS;MBS:RE4; 

meet —tNeeEYS = fF Zi="" THEN LSOELSE TO=VAL (ZI) 

Boe RING gEPt+sl, Z£Zis: FOR =i 7G 100 : NEXT 

meee [8-1 GR IO+KS THEN PRINT EP, EBLSs: PRINT JEP, "2X 
NOT A VALTD OPTION eeek™>: FOR T=1 TG 1000 : NEXT : FRINT 
wEP. ELS; 3: GOTQ 149 

160 RETURN 

Deepen to.tee “st Ss" or "sRes" "sRESs " 

Beans MESs RSS 5 

Pere S—INKEYS = IF 2Z@="" THEN 200ELSE PRINT g@FP+13,2Z0;:: FCR 
T=i TQ 100 : NEXT 

2iQ [TF ZO 5L$ AND ZQ<>CHRS(ASCILS)—-Z2) AND ZOE“ RS AND 
Berens (ASC CRS) —SlbewHEN 1976 

220 RETURN 

ScO REM 

Beets + PRINT TAB(1S) "PROBABILITY DISTRIBUTIONS 
SaQerPRDMGsPRINT: I2=0 

Seer iNi "MENU: " 

340 PRINT" (1) Normal Distribution 

geo PRINT" (2) Wilcoxson Signed-Ranks Distribution" 

moe PRINT ™<¢3) Mann-Whitney Distributison” 

S7O PRINT “(4) Kologmorov Distribution" 

280 PRINT "(3 Siete t osu. bition” 

sage RINT °fS) biliitor’s Normaiity Test Distribution" 
600 PRINT "¢7) Lilliefor*s Exponentiality Test 
Weocrihuticn” 

SeeerRINT "(8S) Chi-Square Distribution” 

620 EP=8976 : KS=8 : GOSUB 139 

6350 ON 10 GOSUB 2000, S000, 5000, 6500, 7000, S000, 9000. 16000 
640 IF FLAG=1 THEN FLAG=0:FOX=1-FOX: OOX=1-00Xx 

Seemermien 5 PRINT 3: PRINT "“Frob. ( X €="sPX3s") = “3s : PRINT 
USING "4+. #4#4#";FC0X 

Geers y "Prob. ( X & "2FX3") = "3 : PRINT 





USING" #. ###H#" 5 QOX 


SoUeeme td we7oO, Run again or Guit?7"s: RP=978 : Ls="r" : 
R$S="g" =: GOSUB 190 
Seema s tees if 7O="Fr" OR 2Z0="R" GOTO SS0ELSE RETURN 


690 GOTO 30 


1000 
1010 


Sesser tne me tio) “INVERSE PROBABILITY DISTRISUTIONS 
PRINT #1925 “MENU: *: PRINT” (1) Inverse Normal 


me2ZO FPreeaNT” (2) Inverse Wilcexson-Signed Ranks" 
£P7S0O PRINT ° (3) Inverse Mann-Whitney Distribution" 
1049 PRINT "(¢4) Inverse Kolmogorov Distribution" 
1050 PRINT "(S) Inverse Smirnov Distribution” 

1960 PRINT "(8) Inverse Lilliefor Normality Test 
Distribution" 

mor PRINT “¢7/) Inverse Lilliefor Exnonential Test 
Beste ibution" 

1080 PRINT "(8) Inverse Chi-Scuare Distribution" 
19090 KP=996:4S=8:GOSUB 130 


11G0 
tee} 


1Z0G0, 14000, 15000, 


PZ 
oO 


PRINT 


1140 
ale 
EL Se 


T2=0 
SNPS SeSUa 


14000, 17000, 18000, 19000, 19220 


mares etrmir) «For Fs) = Prob. ¢€ A <= x } ="3FR 
IF XQ=9ES6 THEN PRINT “x = aunfinity”" ELSE IF X0U=9E-S3é6 


a = Imtinity” 

TF XO=SES6 OR XQ=SE-36 THEN GOTG 1140 
Biles =) se 2 IF ABS (XO) 21000 THEN PRINT 
IF ABS (XO) 7100 FPRINT USING"##4#.#"°:40 ELSE 


ABS (X@)>10 PRINT USING" ###.8#":%O ELSE PRINT 
USING" #4). #434": XO 


Meer INT aS60,° Run Again or Cult"s: LSe"r" s RS=H="q" : 
Pre=o75 : GOSUE 190 

meee i> 7O="R" GR 2O="r" GOTO L1GOELSE RETURN 

2000 FEM * £ * NORMAL DISTRIBUTION « & 

2010 IF f2=1 GOTO 2o7Aa 

eee) CLS:PRINT 220,"NORMAL DISTRIBUTION" 

2030 PRINT: PRINT: INFUT"Enter Mean": MU 

“O40 INPFUT"Enter Variance": VAR 

2050 IF VAR<=0 PRINT"** ERROR xx ":VAR"IS NOT & VALID 
VARIANCE": GOTO 2049 

2060 INFPUT"Enter Normal Statistic teu es Le feo Gren 
PX=NX 

2070 SD=SER (VAR) 

=UGOQ XN=(NX-MU) /SD:AX=ABS (XN) 


ee? AZ—. STBSS2S5AP ( (C-AXFAX) /2) 

2100 NT=1/ (1+. 23516419#4X) 

mt) 

Be Ale 6). 519581 SENT) —(. SS6S6S8SNT02)4+01.781478E8NTOS)-(1.821 


SSGANTO4)94+¢1. 250274*4NTC5) ) 


IF XN<O GOTO 27150 
FOX=1-@0x 

SETURN 

FOX=00% 

Q@OX=1-FOXx 





maf 4) 
SOND0 
Soi10 
ee20 
3030 
SO4D 


Bee PoSitTiIve”: GOTO 


SO50 
S060 
MUST 
S070 
SOS3O 
2080 
SO9O 
eid 
SLO 
Bet 
S24) 


a0 


eaeit) 


2140 
eo) 
2150 
S70 
3180 
2iFo 
=meo0) 
eet 
ae 

Seo 
S240 
ee 0) 
2260 
Se7 0) 
eet) 
Sec0 
2300 


at 

Sa O 
woe’ 
eet) 
SS4O 
S50 
3360 
Ba/7 9 
2 38O 
2OFD 
3400 
S41 
Z420 
S430 
S440 
2450 


RETURN 


oie 2 eee EE CXSON SlsNED-RANEKS DISTRIBUTIC 
CLS: PRINT#iO, "WILCOXSON SIGNED-RAEANKS DISTRIBU 


PRINT: PRINT 
reot “GMfeER wileoxcon Statistic” 
fee oO REM PRINT “#« #« 2 ERRCR 


3030 


ViFX=V 


i Fe | STATI 


INPUT "ENTER Sample Size"3;IN 

tae = ft THEN FRING “siete SZ ERROR £°e & £SAMP 
Be 21 °s6O0TG Sos0 

IF V>INS& CIN+1) /2 THEN FRINT " & kK ERROR x & 
TF V>INs CIN+1) 74 THEN V=IN# CIN+1) /2-V-1:FLAG= 
IF V>27 OR CINGS AND VIN) THEN GOCSUB 2620 EL 
RETURN 

Peet VI —YeGeoun SISO ELSE GOTO S1s0 

RETURN 

V=EINT (CY) 2: GOSUBR 2150: TFOX=FOX: V=V+1:GOSUB 


See oxt A) / 27 80xk=1-FOxX 

RETURN 

fees Ce tite) 0s 12 64) S03 11 (5) =02111(6)=0 
C=1 

N=INT (.5% (-14+S50R(61+8xY))) 

FOR I=i TON 

Ti¢1)=0 

Ti@2)=11 Ci) +1 

Peete. GOEG S360 

peed = 00 

meee) —fhic2)+1 

Peete s GOTO 55560 

wince) —11¢2) 

feet) —iit sire 

IF [<4 GOTG 3360 

Ti¢4)=T1i (3) 

Tics =LTic4)+i 

feet SEITE S360 

Ti¢S)=11¢4) 

eco) — tito) +1 

meet. 6 GOTO 3360 

Tifé6jp=11 (3) 

er.o)—1iits)+1 

eee eet ee yi 1S) 411043411 (S)+11 6S) 
Moroey HEM GOTO 3290 

Peete | CGOC 3S580,35570, 53260, 3540, 3540, 2550 


TF [71 GOTO 34190 
PiL¢1)=IN:s GOTO S580 
SOR TL=1 TO {[-1 


s=9:J=I-IL 
EF J-i<=0 THEN GOTO 
Bem e—1 TO J-1 


O=StTi tk) 


54719 


N kK & 
TION” 


Sele Must 


LE SIZE 
"2f0TO 
SE GOSUE 





S440 
247/79 
S480 
S490 
ea 0O 
Soild 
Seas 
SOUO 
SoA 
Soo 
SueO 
m7 () 
SUBO 
ou70 
S600 
S619 
2620 
SOS0 
3640 
SO6S0 
3560) 
Boy.) 
SéB0 
3596 
3720 
m7 iG 


—=—=-97 =, 
-, / -t 3 
en 


e7s0 


NEXT K 


See lowe Cet £CTi +2) se 


foes ee PHoN GOTO 
mcr i—iLjy=in 


CN 


J GOTG 


NEXT IL 


fie y— PN: GOTO 


ce 
Le 
IF 
a 
Le 
ilies 


T1i¢6)<1N 
PVASY< IN 
Ti¢t4)< IN 
Petey 
PCL ret N 
it). iM 


aes 
as a 


70, 3560, 3550, 5540, 


GOTO 
GOTO 
GOTO 
GOTO 
SOTO 
GOTO 


Sad 


SURO 


SSB @) 


Soe 


se 9D 
3250 


rena Cae heroes 
a ae 


reas 
<3 tI 


NEXT I 

BOA—-C/ 20 IN: GOA=1-FOX 

RETURN 

REM & X * WILCOXSON SIGNED RANK 
MU=INKCIN+1) /4 

VAR= (2N+1) #4 CIN+1)N/24 

NX=V+0.35 

SuSsue 2070 
L4=(SNC2+3N—-1) / C1ON# CIN+1) C2N+1)) 
W= (NX-MU) /SOQR (VAR) SF= (EXP C-O.S4WELS) > 
FS=F*e (WE S-SEW) 

FOX=FOA+L4% FS 

FPSBS=-F Xx (WOS-1LOAWC S41 540) 

a et ey ee WS tt OSAWT S—1OSsW) 


APPREZ DISTISuTICN’* 


PJ Ser(es2oo1S) 


Ll S=4% (SNC4+SN03-SN4+1) / (S58 CINE CIN+1) X(2SN4173 302) 


ee POX=FCAtLSS FE+. SAF7% LACS 

fea OOA=1—-FOX 

3760 RETURN 

S000 REM *« XX * HOANN-WHITNEY DISTRIBUTION x & 

meee CL StPRINTwiO, “MANN-WHITNEY DISTRIBUTION" 

mee FRINT: PRINT 

Seer INPUT “ENTER Mann-Whitney Statistic’ :U:sPx=U 

meee tr UY =< O THEN PRINT “* * * ERROR *€ * © «STATISTIC 
ieee! ££ FOSITIVE *:GOTG T0300 

wOS0 INPUT "ENTER Smaller Sample Size";N 

eet N= 1 THEN PRINT "2° %. 4. ERROR *£°4 & SAMPLE SIZ 
eet BE > 1": GOTO SOs0 

Sey GNPUT “ENTER Larger Sample Size”":M 

eee te M =. N THEN PRINT " & ® kK ERROR * * * YOU ENTERED 
Pee cUeLLER SAMPLE SIZE LAST ":sGOTO 3050 

9090 FLAG=0 

3100 IF USNEM/2-.5 THEN U=N#M—-U-1:FLAG=1 

S110 IN=N+M: V=U+N4 (N41) /2: I=N 

feet UNM THEN PRINT “ & *#°ERROR * * THE STATISTIC 
Sees | THe FRODUCT OF THE SAMPLE SIZES":GOTO 3030 
seme ff 727 OR MSS THEN GOSUB S650 ELSE GOSUE S150 

3140 RETURN 


ero0 


tees P—V hem COSUR Si97O ELSE GOTO S170 


94 





3150 RETURN 


31i7O 
Sleo 
9180 
ai9O 
S200 
we 1 
Sec 
Sa 
w240 
S250 
2G 
west} 
Soi) 
we GO 
ot 
cy 10) 
Soe 
S500 
S540 
Soo) 
350 
Serge. 
S230 
a7 
S400 
23410 
3420 
S450 
o440) 
3450 
Sooo 
54709 
5480 
5490 
3500 
apes, 
S20) 
5530 
SOA 
ott 
BwOU 
sires) 
S380 
ao Fi 
S500) 
3610 
mau 
DOs 
[O40 
Soot! 
3550 


V=INT (V) :GOSUB 


FOX= (FOX+TFOX) /2: GO*A=1-FOX 


RETURN 


Ti (2) =02 1163) =027T104) =0:7T1 0S) =O: 7T1C6)= 


C=0 

Ti¢id=0 

Pec y=ficiyst 

Ne i<2weOTG 3370 
tail) 
eve — tit 27d 

fa i-S GETOG S380 
el (3S) =11 (2) 

Pees (S113) +1 

ee L24 GOTO 3360 
Ti ¢4)=11(3) 
Ti¢4)=11¢(4)+1 
Paes COP eo 5eO 
Tits) =114) 
Pitad=li (a) +2 
Poa GOP Ss le0 
I1¢6)=f1¢5) 
Pavov—ilic) +1 


IF S?V THEN 
C=C+1:0N 
IF I31 GOTO 3436 
I1¢1})=IN:GOTO 
FOR IL=1 TO I-1 
S=O: J=I-IL 


SOTO 


Be 3-1. =0 THEN GOTG 


FOR K=1 TO J-i 
S=S+11(K) 
NEXT KE 


SaortIitiy itd) +c! 


beeos/Y THEN GCOTG 
Pee 1 — TL) IN 


Eiger ocu.e Ssoa70, S280. sa7 


w410 


aE600 


et 
Soe 


NEXT IL 

Tift =IN: GOTO S600 
foo?) met GOTO SS70 
Poe aN SET S340 
tet boy Pee GTO. S510 
Ie LitS)< ieee TO S280 
Hata IM GOTO S250 
Mette IM GOTC 3220 


S=1 


pene) 1 Note S=SeeGIN—-k) / (N-E) : NEAT 
FOX=C/S: QOX=1-F OX 


RETURN 
REM 


 * * MANN-WHITNE 
/ 


MUHS=N/22VAR=N#EONtM+LG 


7 


2490 


if 


FR 


ie 


PPROXIMATE Di 
22S 


TSESGTHO SSE00., 5590, 5580 


I Pies) ve 


C4 


R 


= 


QO 


TR 


sess iH OA—- OX: V—V+1> GOSUE 


"od 


fe? 


7 eit 
ie dee” 


af OS 


eas «aaa 





Sev ceotn 257) 

S680 L4=(MC2+NO2+ITN4+M4N) / (204AN# CINAT) ) 

Stele A®, 

1 S6=(2¥ (ML44+N0D4)4+44N4 (MI S4NES) +68 (MOS) KONE SD) +44 (MD S4+NOS) +7 
¥NS CMAN) MC S4+NDS+2aN-TND SC 21OX (MES) EINES AC CINE) OD) 
peed We PASM 7 SER (VAR) sPH (EAP (-O. 502) )/S0R(6. 28318) 
et) fS—F £ (WE S-S) 

wee) Po=—-F# CWE S—-1OCS+15) 

epee 7 ——F 24 7-27 £S+1 0503-105) 

o/40 FOX=FOX4+L4& FS 

eaea FOX=FOX+LGs FS 

S750 FOX=FOX+.5% F7k L4AC2 

2770 RETURN 

S000 REM #*£ &* ® FEOLMOGROV DISTRIBUTION x x 

Seto ClLS:PRINT310, “KOLMOGOROV DISTRIBUTION" 

eucO PRINT: PRINT 

eee INPUT “EMER Skolmcgorov Statistic";D 

60460 IF DLO GR D?1i THEN PRINT "“X * ERROR & £¥":GOTO 6020 
6950 INFUT “ENTER Sample Size"s5N 

£060 IF Neil THEN PRINT " # ERROR #€ £€":GOTO 6tho 
6070 S=D4(0. 12+SCR(NI40.11/SGR0N)) 

65086 GOSUB 6110 

5090 PX=D 

6100 RETURN 

epee OOX—C. LSS) CCSL62) 

6120 IF GOX<0O THEN @0Oxk=0 

6130 TF Q0O471 THEN GOxX=1 

6140 FOx=1-80xX 

5150 RETURN 

7000 REM & & * SMIRNGY DISTRIBUTICN x # 

Pee CLS:FRINTS10, “SMIRNGY DISTRIBSUTICN" 

mee PRINT: PRINT 

fee INPUT “ENTER Smirnev Statistic":D 

eee tr 560 OR OFl THEN PRINT “*£ * ERROR *« *":GOTO 7O30 
7OS0 INPUT "ENTER Smaller Sample Size":N 

eee ete Nel THEN PRIN? “* & ERROR *&® * SAMPLE SIZE MUST 
oe 2"°S50TO FSO 

Bee eINPUT “ENTER Larger Sample Size": 

feeoe i> MN TREN PRINT “*' * ERROR # x YGU ENTERED THE 
Seg Serle Size LAST":GOTO 7050 

Peo Te Ne=.i THEN S=(D—-1/(2)14(.i12+SGR(N)+.11/SGR(M)) 
LSE ITF (M/NO=INT{M/NM) THEN 

S=D# SGRCN/ (N+M))4+2/ (3% SQR (MM) 3} 

S=D* SOQR(N/ (NM) )+2/¢(5% SORUM)) 


ea 


b 


{ 


on 


rm 


pare GOS 
7110 PX=D 
7120 RETURN 
Seeeere!s 45% *@ j\LILLIFOR’S TEST FCR NORMALITY «x & 
8010 J=I 

eee ease RIM] “SELECT CPTION'’:PRINT:FRINT 

eeew PRINT “¢1) MEAN and VARIANCE are unknown": FPRINT 
ee) PRind “ (2) a is known, VARIANCE is unknown" 


7& 





8050 KFP=G96:K5=2:60SUB 1406 

SeoeeeeogemimiaiG, "LILLIFOR’S TEST FOR NORMALITY" 

8070 PRINT: PRINT 

S080 INPUT "ENTER SkKolmogorov Statistic"s;D 

Seow t- 2.0 Ge Dri THEN PRINT “4 % ERRGR &«£ %":GOTO scso 
8100 INPUT "ENTER Sample Size";N 

8110 IF N<i THEN PRINT "« * ERROR x # SAMPLE SIZE MUST RE 
Bier tIve’=:GOTG S100 

8120 DENOM=1.0:PA=D 

8130 D=Px%/DENOM 

Si4o IF (G=1 THEN S=D% (-—.01+S0R (N)4+0.95/SGRiN}) ELSE 
S=D#¥ (0. O35+S0R (NI+.2983/S0R (N)) 


8is0 TDENCGM=DENOM 

eeol GOSUB 451190 

See iF IGFi THEN GOTO 82560 

Sree if FOX <= .85 THEN DEN@M=1-.24¢ FOX:GOTO S250 

meee i= FP@A <= .90 THEN DENOM=1.357- .464 FOA:COTO 8230 
Beet FOR <= .FS TREN DENOM = 1.359 —- .644% FOX : GOTO 
S25) 

oes x a= 6975 THEN DEMGM = 1.795 -— 1.12% FOX =: EATO 
Seo) 


S220 DENOM = 2.1335-1.467xk FOX 

S230 IF SRBS(TDENOGM-DENOM) > .801 THEN GOTO 8130 
S240 TO=J 

S230 RETURN 

S250 IF FOX <= .9 THEN DENOM=1.85-.8% FOX :-GOTO S230 
=e7o I[F FOX Peo SeeM WeNGM =1.47 — ~.4x% FORK =: SO 
~» o 


5 = y 
Seo [fF FOX < See weMeiy = toa o — ~ek FOX 6G Bes2e5u 
2 


=. Saou 
geo! DENOM = 1.07 SOTG 8230 

Seite € £94 jLILtLIFOR’S TEST FOR EXPONENTIALITY «* ™® 
Perec Ss PRINGeiO, “LILLIFOR’S TEST FOR EXFONENTIALITY" 

pom PRINT: PRINT 

Seee tTNEUT ®ENTER Skoimecgorov Statistic”";D 
Geert D220 OR D>t THEN PRINT “* 0% ERROR & 
wees INPUT “ENTER Sample Size":N 

Beer tt N<i TREN PRINT “2 x ERROR #£ * SAMPLE SIZE MUST SE 
Beet FiVe™:6GOTCG 8100 

Beet © + 1.3442/ (SOR (N)+.256+.5/S0R(N))+.2/N THEN FR=1 : 
G07 9140 

elaee PR=.3:50SUB 19090 

e799 IF DeXG THEN P1i=PR:FP2=1:GOTO 91146 

F100 Pi=0:FP2=PR 

fee PR=(P1i+P2)/2 

eee COSUB 19090 

i130 IF ABS (D-X0)<.0001 THEN GOTO 9160 

Sta) IF D-X@ THEN Pi=PR ELSE P2=FR 

meee GOTO Srro 

150 FOX=PR 

Piya OO%=1-FOX 

mee PX = YD 

geet! RETURN 





10010 
1o02zD 
1OQS30 
Pex 

10040 


Mwot BE > 


10050 
19060 
19072 
ENTER 
10086 
19N90 
10100 
ae 
eT O 
1Q130 
fOL40 
Pero) 
10166 
gL 7 O 
10180 
TH1LSa 
Loto 
1931 
10220 
1O2 20 
LOoz4a 
19250 


102 
1o270 
20280 
LOe9O 
16500 
feu 1) 


oO) 


10080 


Nie weo. sent —-SeUARE DIST ee 
ie oie 3192, ""3:INPUT"Enter Chi-So ee Seat suis 
IF CX<=0 PRINT"X% ERROR 4% GHESSGuane, STATIST 


wre TO 1eGs0 

TF [2=0 THEN PX=CX 

INPUT"Enter Decrees of Freedom"; DF 

IF DF<=0 OR DFA+INT(DF) THEN FRINT"#%4 ERROR 4 


AN INTEGER + O°:50TO 10060 
teem so Cle fos 70 

DP=1 

per —fitooe) 10 2 STEP —2 
DP=DP 

NEXT I 


MeE—i ni tne (CDF +1) /2))#E 
it INT (DF/2) =DF/2 GOTO 10170 
Beas As Se LoL SFG) 

COP Vi eled 

PC=1 


D=D+2 

MC=NCX CxX/D 
Pee OUI Oe T 
EC-Ee She 
CG 1021 
FOX=PCANPX 
Ho,=1—Fex 
RETURN 

in (We a C1 a) 
MU=9 2: VAR=1 : 
EESUR 2070 


ee 


(1-(2/ (94DF}))) /SOR(2/ 


ve 


{Swape j.5 
. - °a* tof 2 


a od, 


1O220 RETURN 

Pomeo REM <<) & = INVERSE NORMAL DISTRIBUTION $$& & 

LS010 ff T2=1 GOTO isiac 

Meee! CCS =: PRINT 320, “INVERSE NORMAL” 

Peet PRINT = PRINT =: INFUT"Enter Mean";MO 

15040 INFUT "Enter Variance"; VR 

5050 IF VRL=0O PRINT" *X ERROR &x ere. ko ST ee VAL TD 
YARIANCE":GOTG 13040 

Meeoe INEUTY Enter Frobability (1.8. Fix} = Frob. ¢ X <= 
me) eR 

13070 IF PR=0 THEN XOQ=9E-S36 ELSE IF PR=1 XO=964+36 

12060 IF FR=0 GF PR=1 RETURN 

15090 IF FR<O OR PR?1 THEN PRINT"x* ERRCR £2 afta oe ey 
A VALTD PROBABILITY’: GOTG i23060 

mee I PRLS THEN RL=SOR (LOG (1/(41-PR)0T2))) Bele 
RL=SOR(LOG( i/ (PRIS) )3 

13110 

AG=RL—f (2.3515517+(. SO28SS4RL)+¢(. C1 OS2SSRKRLI2Z) 3 /(14+(1.43527S9EGkR 





ks 


X= 


1892S9HRLI2) +6 
13120 IF FPR>= 
MO— (SOR (V 


SeG4+REES))) 
XOS=MO+ (SOR CVR) XO) 


Seen 
«Sane! 
RD XO) 


i3120 RETURN 


14000 


Pee CES. eRING 


RES Seege'ses INVERSE 


Bees? lB) TON” 


14020 
14030 
14040 


poem sP RINT 
INPUT "ENTER Sample Size"; IN 
iP iie< 1 THEN FPRONT “ax X 


feot BE POSITIVE":GOTG 14030 


14450 
1404609 
14aq7a 
14071 
14072 
14073 
14074 
14075 
149785 
14980 
14090 
14100 
t4119 
Hees? it) 
pate) 


MO=-INS CIN+1)/4-.5 

VR= CA2N+1) 0 IN+1L)N/24 

SESUB 12060 

V=INT (xX@) 

GOSUB 3090 

TFOX=FOX 

V=V+1 

SESseweuss 

XGQ=INT (XO) +(CDBL CPR) -TFOX) / (FOR- 
MOSAG-INT (XQ) 

TP MO... 25 THEN AUSINT (AG) s RETURN 
Pee soe TeeN KAO-INT(xAO)+.5 ELSE 
RETURN 

REM. ee Ke INVERSE 
GESsrr Lt " INVERSE MANN 
a se NY 


WHITNE: ¥ 


aed 


ERRORS 2 x 


DIS 


ee 


WILCOXSON SIGNED-RANKS x # 
"INVERSE WILCOXSON SIGNED-RANKS 


SAMPLE 


XO= CA 
* 
im bee oN 


13020 INPUT “ENTER Smaller Sample Size":N 

iSC40 TF Ni THEN FRINT "XX x & ERROR * * # SAMPLE SIZES 
Seeeocii Ive’: GOTO 15030 

meeoo INPUT “Enter Larger Sample Size":M™ 

feos IF M1 THEN PRINT "* * X ERROR * =x x SAMPLE SIZE 
peweesSsirive’:sSOTG 15020 

15G70 MO=N4M/2 

130680 VA=NaMx (N4+M41)/12 

meer GOSUB 12060 

15091 U=INT CX) 

meee OOSUB S130 

15093 TFOX=FOX 

15094 U=U+i 

meecS CGOSUB S150 

mes? AG=INT (AG) +( CDS (CPR) -CDBL(TFOA)) / (FOX—-CDBLCTFOX) } 
15100 MC=*YQ-INT (4B) 

i3110 TF MO@<.26 THEN XOQ=INT¢(X@) : RETURN 

15120 IF MOz.7=S THEN 4OQ=INT(XO)+.5 ELSE XO=INT (XO) +1 
1351230 RETURN 

15000 REM «k * &¥ INVERSE KOLMOGOROV « # 

15010 CLS : PRINT 920, "INVERSE KOLMCGORGY" 

meee PRINT: PRINT 

Bees? INPUT" ENTER Sample Size":N 

16040 IF N<1l THEN PRINT "x * ERROR & £":GOTO 14030 
Meese INPUT enter FProSability (i.e. F(x) = Frob. ¢ X 


SIZE 


ee hom 


C 
Ui 
“4 


wa 4 >, 
mus 





eo eR 

16060 IF FPR=1 THEN X@= 0.0 ELSE IF FPR=0 40=5E+35 
164070 IF FR=90 OR FPR=1 RETURN 

14080 IF PRAO OR PR?1 THEN FRINT“&** ERROR #2 meseiaee koe oC), P 
A VALID PROBABILITY’: GO0TQ 16020 

Peeero S—-SGR(LOG(I-PR) /LOG(.135) >) 

16160 D=S/(.12+S0QR(Nji+.11/SQR UN) }) 

ié1i0 xX@=D 

16129 RETURN 

1790090 REM * kK kK INVERSE SMIRNGYV & & 

meet CLS 7s FPRINT 29320, "INVERSE SMIRNOV" 

17920 PRINT: PRINT 

17030 INPUT"ENTER Smaller Sample Size":N 

176490 IF N«<il THEN PRINT "* xX ERROR € *4°:GOTO 17920 
17050 INPUT"ENTER Larger Samele S122"; 

17950 ITF°OMAN THEN FRINT “*% ** ERROR * #£"°:GOTO 17930 


mee? INPUT Enter Prebabitity ‘i.e. Fix) = Prob. ¢( % 4= x 
) = "SFR 

meood IF PR-1 THEN *0= 0.0 ELSE IF PR=0 THEN AQ=9E+26 
170909 IF PR=0 OR PR=1 RETURN 

17100 IF FPR<O OR PRi THEN PRINT" &% ERROR 22% fers 1 Sao 


A VALID PRCBABILITY"':GOTOQ 17060 
17110 S=SOR(LOG(1-FR)/LOG(.135)> 
$7120 IF N<=0.14M THEN D=S/(.12+SOQR (N)+.11/S0R(N))4+1/ (24M) 
ELSE IF (M/N)=INTCM/N? THEN 
={(S-2/ (3k SG@R(M))) /SOR CNEM/ (N+M) ) 
Fi SE D=(S-—2/ (3% SOR(M) ))/SGR(UNSM/ (N+M) ) 
feito) AB=D 
17140 RETURN 
18000 REM * * * INVERSE LILLIEFOR’S NOGRMALSLITY x *¥ 
fee J—IGSCLS:PRINYT “SELECT OPTION": PRINT: FRINT 
18020 PRINT "(14> MESN and VARIANCE are unknown’ PRINT 


meso PRINT " (2) MESN is known, VARIANCE is unknown" 
18040 KRP=895:K5=2:GOSUB 140 

PoeeriseosreiMNs wis, “INVERSE LILLIEFOR’S NORMALITY TEST 
DISTRIBUTION” 


19060 PRINT:PRINT 
19070 INFPUT"ENTER Sample Size":N 


mee Te N<i THEN PRINT “sk # ERRECR 4&4) #":GOTOQ 16020 

foe Ne emter Probabitity (1.8. Fix) = Prob. ( K <= 
m= 86h R 

18100 IF FR=0 THEN X@=98-36 ELSE If PR=t X8=9E+=26 

18ii0 IF FPR=0 OR PR=1 THEN IO0Q=J: RETURN 

poe 1 PR<O OR PRE1 THEN FPRINT"*£# ERROR xx his to er 
Pe YALID PRORABILITY”"’:GOTO 130606 

meeeo S=SUR (LOG i (i—-PR)) /LOG(.135)) 

18i40 IF IG=1 THEN D=S/ (-.G1+SQR(N)+.85/SOQR(N)} ELSE 

D=S/ (. 055+. 283/SQR (N) +S@QR (N) ) 

Mees? fF 1O>1 THEN GOTO 16240 

meee I> PR <= .35 THEN DENOM=1-—.24® FOCX:GOTO 18210 

fae im PR == (90 THEN DENOM=1.357-— .66*PR:SOTO 182106 
fem if FR 2= 7-975 TREN DENOM = 1.3397 - 0.64#FPR +: GOTO 


os, os 
a3 i J 





iG 
Poeeceee — ee S PEN DENOM = 1.793 — 1.12sPR =: GOTO 


P3200 DENOM = 2. 1S35-1.447%kPR 

18210 XG=DxDENOM 

tWee20 10 = J 

18230 RETURN 

Peete ee == 29 THEN DENMOM = 1.85 - .S4PR =: GOTO 18210 
Pees 8 PR «= 985 THEN DENGM = 1.47 — .44PR : GOTO 18210 
Meeoo ie PR <= .975 TREN CENOM = 1.85 — .84FR : GOTO 18210 
hoor DENGM = 1.07:G0TO 18210 

Pees REM oF i wae INVERSE LILLIEFOR EAPONENTIAL TEST xX & 
epee eioseeiNt 2iS, “INVERSE LILLIibFOR’S EXPONENTIAL TEST 


S ff Ae 


é 


| 


19020 PRINT: PRINT 

19020 INPUT"ENTER Sample Size":N 

19940 IF Ne<i THEN PRINT "* * ERROR * *":GOTO 146020 

Peo? INPUT Enter Fropability (i.e. Fx) = Frab. ¢( X <= x 


} = *-ER 

¢ 4 e * 

19050 IF PR=0 THEN X@=9E-36 ELSE IF FR=1 XO=SE+2& 
19070 IF FPR=0 OR PR=1 RETURN 


Booed IF PRO OCF PRI THEN PRINT"4£% ERROR £3 “eins o> Gene 
Seeveei® PROBABILITY" :&6070 19050 
fees? DENOM=SCR (MN) +. 25+.5/S0R (N) 


Pe Oe ee ee eee AC! = So. 1 2SeeR — 3.7508 3; SOTO 17199 
Beem te Pa —> .7C THEN A@ = 2.0264PR — .8367 : SCTOQ 19190 
meee te eri => .80 THEN AQ = 1.1242FPR -— .0249 : GOTO 19199 
meee tee eR => .7 TREN XO = ./74S8PR + .2799 : GOTO 19190 


meee IF PS => .35 THEN 4@ = .Ss46¢PR + .4187 =: GOTO 197190 
Pee iP PR => .3S THEN AQ = .46544PR + .4583 : GOTO i9199 
Meee if PR =F .2 THEN xX@ = .4852PR + .4525 3: GOTO 19190 
Peewee te FPR => .1 TREN A@ = .S8S72PR + .4521 3: GOTO 197199 


197180 XO = .8464PR + 1.4062 
19190 xXQ=KE/DENCM 

Meeee AD = KO + .2/N 

feet) RETURN 
Peee) REM 
fees IF Ie 


Bve0o0 CL 


Pm oO Ge. 2070 

Beta e BtHVe RSE ECHI-SOQUARE” 

2U0O10 PRINT : INFUT"Enter Degrees of Freecom": DF 

meee TF DF<=0 OR DFAFINT (DF) THEN FPRINT"&£3 ERROR £2 

ENTER AN INTEGER > 0":GOTO 20010 

29050 INPUT “Enter Probability teen ok? = FProo. 1 A 1S x 
me +P R 

20040 IF PR=0 THEN XQ=0 : RETURN 

peo TF PR=1 THEN XQ=SES6 =: RETURN 

eeeoo IF PRO GR PR+1 THEN FRINT'**s ERROR x& eee eo 
A VALID PROBABILITY" :GOTO 20030 

2V070 MO=0 : VR=1 : I[2=1 : PZ=FR 

eee IF DOFri GOTO 2o1s0 

Sets) PR=.oPLL—PZ) 

29100 GOSUB 13100 





S[O110 XO=xXQC2 3: PR=P2Z 

20120 RETURN 

Seee. iP pees GOTO 20160 

20140 X@=-2x LOG(1-F2) 

eego0 RETURN 

eetow GOSUB 15100 

20170 DQ=xXQ 

20180 IF DF?>(2+INT (4xABS(DQ))) GOTO 20280 

Seno — 1 POOR — 2 22S/7S6e/DF — .01513904/DFL2 

meee eo —ATEIDP Elo 6 DO Kk (.4715941 + .02607083/DF —- 

moms TSGOOs7 /DFL 2) ) 

eeeeo *9=KS + t1/DF * DOCS X (€.00013480278 + .01128186/DF + 
meer? o/IS/DFLE 2) } 

ree eee Ao CUE Pil SG DELS ek (-. COESSS069 -— .011535751/DF 
meee oe so / OF C2) ) 

Peewee ASAP + (DFi-=2 * DELS * (.00S312558 + .005169654/DF —- 
Bee C SOS 5G/DF £2) ) 

Bees AG-AF + (DPE-2.5 * DELS * (-—. 0008426812 + 

meee SOOL/DE + .0O01060438/DFL2) ) 

Sees e As —A7 = (DFi-o *£ POLE pe 6. 00009780499 -— .0014t0117/0F 
Pee aoa eo/ DP LT) ) 

2260 XO=DFRASES 

20270 RETURN 

eae ear —t— 2/194) + (4e0D004 + 16eDOCS —-2B) / C1 2SLSaADFCS) 
2O=SZ9G X9=K9 + (BkKDGLS + F2OKDOL4 + SILVsS#*DOCS + 
meet? / i 2 eS GSaXRPFLIS) 

eee S—bO/ 5 + i -DOLS + SADQ) / (1 S220F) 

20310 *B8=x%8 — (S*DOQLS + 404D0C3 + 45kD0) / (SSSR2RDF C23 

PUSS? XAS=xS + (3014DOQL7 -— 15174DQCLT —S27E69xuDQ05 —- 
fFSAIGKDO)) / 0787 S200KDF CS) 

M@oo0 AS= SOR(2/DF) * AS 

pees) 2O=UF & ((X84+4A9) £35) 

oes RETURN 
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